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HAUPTAUFSÄTZE 


Geschwindigkeits- und 'lemperaturverteilung hinter einem 
Gitter bei turbulenter Strömung. 


Von R. Gran Olsson in Trondheim. 


1. Einleitung. Die Untersuchungen zur ausgebildeten Turbulenz sind in neuerer Zeit 
durch Arbeiten gefördert worden, die auf einen Ansatz von L. Prandtl über die scheinbare 
Schubspannung zurückgehen). Insbesondere hat sich der Ansatz im Falle der sogenannten 
'reien Turbulenz, bei der also keine Wandreibung auftritt, als sehr fruchtbar erwiesen. Da 
(die Ableitung des Ansatzes an verschiedenen Stellen in der Literatur gegeben wurde, soll 
hier nieht näher darauf eingegangen werden. Für eine Strömung in der .r „Ebene parallel 
zur „Achse mit einem Geschwindigkeitsgefälle quer zur Strömungsriehtung ist die schein- 
bare Schubspannung | 

„ou ou 
wo «die Geschwindigkeitskomponente in der „-Riehtung und o die Dichte der Flüssigkeit 
bezeiehnen. 1 ist der sogenannte Mischungsweg und spielt hier eine ähnliche Rolle wie die 
Weglänge in der kinetischen Gastheorie. 

Vergleicht ınan diesen mit dem für die Sehubspannung bei gewöhnlicher Zähigkeit 

2eitenden Ausdruck 


—MU 
> wird die der Zähigkeit z. entsprechende Austauschgröße 
„ou 
- 


© Größe, die dem Austauch unterliegt, ist der Geschwindigkeitsvektor, der hier dureh seine 
"wiegende Komponente « ersetzt Ist.) 
Im Jahre 1932 veröffentlichte G. I. Taylor einen Ansatz, bei dem die Wirbelstärke der 
'ptbewegung die Größe ist, die dem Austausch unterliegt, während der Ansatz von L. Prandtl 
‚L. Prandt1:ZAMM, Bd.5 (1925), S. 136 oder Verhandl. 2. Intern. Kongr. Tech. Mech. 1926, Zürich u. Leipzig 1927. 
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so gedeutet werden kann, daß der Impuls ausgetauscht wird’). Taylor gibt nicht dir 


O7 
die scheinbare Schubspannung 7 sondern 77 an, nämlich 
Or 


Für T==konstant kann dieser Ausdruck bei positivem \ - 
( 


zu 
B (° u 

\0y 
integriert werden, unterscheidet sich also nur durch den Faktor 5 vom obigen Ausdruck für r. 


Aus dem Ausdruck für die Austauschgröße ® geht hervor, daß der Austausch an Stellen 
der Gesehwindigekeitsmaxima und -minima verschwinden müßte. Dies ist in Wirklichkeit 
nicht der Fall, weshalb L. Prandtl eine Verfeinerung des Ansatzes für e angegeben hat’) 


. . . ( 
Danach hat man den Austausch proportional einem statistischen Mittelwert von y 


über eine 


(du 
kleine Umgebung der betrachteten Stelle, also proportional | (5 > 


gibt den zeitlichen Mittelwert an). Hierfür darf man angenähert schreiben 


zu setzen (der Striel 


- 


wobei die Länge 7 im allgemeinen dem Experiment zu entnehmen ist. Die Austauschgröße 
hat dann «die Form 


- 


(du? 


Mit diesem Ansatz gestaltet sich die Reehnung meist viel schwieriger als mit dem einfacheren 
Ansatz, jedoch ist in dem vorliegenden Fall unter einer bestimmten Annahme über I’ die 
Rechnung sogar leichter durehzuführen als mit dem üblichen Ausdruck für den Austausch. 

Im folgenden soll nun die Gescehwindigkeitsverteilung hinter einem System paralleler, 
äquidistanter Stäbe für 7-0 und 7==0 berechnet werden. Weiter wird die Temperatur- 
verteilung bei geheizten Stäben, ebenfalls für beide Annahmen über den Austausch, angegeben, 
aueh im Fall eines geheizten Stabes. Im Abschnitt IT werden die Meßergebnisse mitgeteilt 
und mit den theoretischen heehnungen verglichen. 


l. Iheoretischer Teil. 


2. Geschwindigkeitsverteilung. Denken wir uns jetzt ein aus zylindrischen, sehr langen 
Stäben in konstanter Entfernung bestehendes Gitter senkrecht zur Strömungsrichtung in einen 
Luftstrom mit gleichmäßiger Geschwindigkeitsverteilung gebracht, so liegt eine zweidimen- 
sionale Strömung vor. Die Aufgabe ist, mit Hilfe der Prandtlschen Ansätze im turbulenten 
Vermisehungszebiet hinter dem Gitter die Geschwindiekeitsverteilune zu bereehnen. Die 
Breite des Vermischungsgebietes nimmt hinter jedem Stab zu, bis sie gleich dem Gitterabstan« 
), wird. Da der Mischungsweg I proportional dieser Breite gesetzt werden kann, so erg) 
sich in genügender Entfernung hinter dem Gitter T= konstant. 


. 


Abb. 1. Festlegung des Koordinatensystems. 


°:,G.I Taylor: Proe. Roy. Soe. A.. Vol. 135 (1932), S. 685 oder auch F. W. Durand, Aerodynamie The 
Vol. III, S. 175, Berlin 1935 (Beitrag von L. Prandt)). 
>», L. Prandtl:a a.0.S8. 138 oder H. Schlicehting: Ing.-Arch. 1 (1930), S. 541. 


4 4 
wobei die Austauschgröße # aus dem Prandtlschen Ansatz übernommen wird, also 
Or 1: ou 
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; dire 4 Die -Achse verlaufe parallel zur Symmetrieebene des Gitters senkrecht zu den Stabachsen, 
‚Achse durch die Achse eines Stabes im Gitter und parallel zur Riehtung der ungestörten 
mung (Abb. 1). Die Geschwindigkeitskomponenten seien « und », die Geschwindigkeit 
ungestörten Strömung U. Hinter dem Gitter sind (in genügender Entfernung) die Ge- 
windigkeitskomponenten: 
u U- v—=Vv,, 
hei u, und v, klein gegen U sind. Man geht nun so vor, daß man in den Bewegungs- 
» schungen nur die Glieder von der niedrigsten Größenordnung berücksichtigt. Unter der 
snahme, daß das Druckglied zu vernachlässigen ist, erhält man eine Differentialgleichung 
ur 10. (Diese Annahme läßt sich nachträglich an Hand der Lösung rechtfertigen.) 
a) Aufstellung und Integration der Differentialgleichung für "=E0. Die 
bowegungsgleichung für die «-Richtung lautet: 
[A (1). 
k für r. Führt man nun U ou, wobei u, < U, U, so erhält man in erster Näherung 
Stellen lie Gleichung 
| ( Or 
lichkeit U, 2 
n hat’). o04 


has Sina Wegen des periodischen Verlaufs von «, liegt es nahe, für u, anzusetzen 


trieh wo «a so bestimmt werden muß, daß die Periode mit dem Gitterabstand 4 übereinstimmt, 

d.h. «a 
/ 


Für ergibt sich dann 


0 
| | ‚| = A’Pa?sin cosay. 
Setzt man nun hY=—=,-, so wird 
4 a 
'acheren 
stauseh. 
ralleler, 
peratur 
regeben, voy 
itgeteilt Setzt man dies in die Differentialgleichung (2) ein, so ergibt sich 


U. 


Da 7 konst. sein soll, so ergibt sich 


in einen 


eidimen- d.h. 
bulenten 
= a Damit Jautet der Ausdruck für die Störungsgeschwindigkeit in erster Näherung 
abstant 
Geschwindigkeit », der ersten Näherung erhält man aus der Kontinuitätsgleiehung 
ou 
| 
Ö u, 1 ] lt | 
= — Ax’’cosaı erhält ma 


H 
2 
r 

S 

li 
4 
. 
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— 


4 


Die Annahme =, „ ist, wie sich zeigen wird, keineswegs so willkürlich, wie es im ers n 


Augenblick den Anschein hat. Mit Hilfe dieser und der für "=-0 geltenden und noch »u 
ermittelnden Lösung werden wir Näherungslösungen für andere Werte I’ aufbauen. \ir 
gehen daher zur Erörterung der entsprechenden Lösung, für !=0 geltend, über. 
b) Integration der Differentialgleichung für Im Falle erlilt 
der Ausdruck für die scheinbare Sehubspannung die Form 
ou om 
. 
0 Q Y Y 
Die Differentialzeleichung lautet 


21° 


wo k 7 ‚ Indem 7, statt 7 gesetzt wurde, da nicht zu erwarten ist, daß der Mischungswog 


bei «den verschiedenen Annahmen über r mit demselben Wert herauskommt. 
Setzt man nun ähnlich wie im Falle 


N A, u. 


und zwar soll F, (0) 1 sein (J, ist also wieder die Amplitude), so erhält man folgende Diffe. 
rentialeleichung für FF, 


Fy)=AkF, 
Setzt man ferner 9 en, wo e eine Konstante ist, so ergibt sich 


F (en)= F,(en)F, (cn), 


wo und (e Jetzt nach » differentiiert sind. Da e willkürlich ist, so wird | 
gesetzt, womit die Gleichung lautet 

Diese Gleichung wird mit 7, multipliziert 
und einmal interriert 
2 
Fi 
| (1— 5) 
g 
wo y=n, für FR =0. Da für F=1 wird, 
dF, 2°\ 3, 
2) 
0 
indem allgemein?) 
| 
dt 
n) 


%, Siehe z. B. E. Jahnke-F. Emde: Funktionentafeln, S. 95, 2. Aufl., Leipzig u. Berlin 19393. 


hr. f. ang u t5 - 
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— 


Im ers on | | 

.0,88623 , | -—)!=135412,  —!=0,927 

en. Wir / / 

sieh 7, = 1,1300°). Da für FR =0, 
erlilt Fi 
dF, 


(Ss Integral |; _ Fa läßt sich z. B. mit Hilfe der Simpsonschen Formel berechnen. 


ist, so wird e: 


/ 


N 
einander entsprechenden Werte F, und 7 bzw. y sind in der folgenden Tabelle I an- 
schrieben 


Tabelle 1. 


F, cos 2a 4 7 F, cos 27 / 
7 
hungsweg - 
1.130 0.250 0,0 0.000 V.AT6 0.105 0.7 0,780 
1.043 0,231 0.121 0.361 0.877 
0.055 0,212 0,2 0,241 0.226 0.050 0,0 0,051 
0,865 0,192 0.3 0,360 0.017 0.98 0.004 
nde 
0.774 0,171 0.4 0.475 0,047 0.010 0.09 
0.680 0,151 0.5 0.586 0,000 0.000 1.00. 1.000 
0.581 0.129 0,6 0,691 
In der letzten Spalte sind die F, entsprechenden Werte für U angeschrieben. Die 
relativen Abweichungen betragen bis zu 20°/,, die absoluten bis 9... 
l,-k 
„3 | Da nun 
U 
22°4580° 
wird 
6 
Die Amplituden müssen in beiden Fällen (für 1’ = Sp und "=0) übereinstimmen, weshalb 
wir z. B. /, dureh ausdrücken können 
| 
Sa’\ı) 2.4,520° 
woraus 
a Um eine brauchbare Näherungslösung für beliebige Werte 7’ zu erhalten, wird folgende 
Überlegung gemacht: Die Störungsgeschwindigkeit u, setzt sich aus der Störungsgesehwindig- 
keit für und der Korrektion, die hinzukommt, wenn zusammen. Es läßt sich also 
Ansatz machen 
/ 
wo gleich «, für ’=0 und die Korrektion. Für "ecosay, 
vonmit wir für @,,, erhalten 
(da). °) Die Werte der Fakultäten sind entnommen: K. Hayashi, Fünfstellige Funktionentafeln, 8. 158, Berlin 1930. 


Ineiner in ungarischer Sprache ersehienenen Arbeit von E.Anderlik (Mat.termeszett. »?,54bis79, 1034) 


rselbe Fall (für ? 0) gereehnet worden. Das graphisch ausgewertete Integral (a) gibt Anderlik zu 1,2824 an, 
end der genauere Wert 1,2936 beträgt. Der Fehler (etwa 1%/,) ist belanglos. Die Arbeit von Anderlik wurde 


‚erf, erst bekannt, nachdem diese Rechnung abgeschlossen war. Dureh die Substitution (1 Fı°) Ia z läßt sieh 
'egral (5) auf ein elliptisches Integral von der Weierstraßschen Normalform bringen. 


| 
_ _ U U U U U U U U U U U U 
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Also ıst 


Diese Formel ist für P=0 und =, „ genau. Für andere Werte von !’ ist sie nur als In! 


polationsformel zu betrachten, sie dürfte aber für O<l’< Se ziemlich genau sein. Da (ie 
bisherigen Messungen allgemein auf einen kleinen Wert 7’ schließen lassen, so dürfte der tat. 
liegen, der von der angegebenen Form) 


sächlicehe Wert von in dem Bereich lo. 


- 


erfaßt wird. In der Abb. 2 sind für "=0 und ’: 5, die Kurven der Geschwindigkeit. 


verteillung aufgetragen. 


005 00 0% 020 050 035 00 045 050 


08 0 00% 00 0% 00 025 090 05 0m 06 0% 
1. A5823] Y 
1.0 
Abb. 2. Geschwindigkeitsverteilung max Abb. 3. Geschwindigkeitsverteilung max 
für ’=0 und!’=,—. für und I!’ 


Die Geschwindigkeit », für 7-0 erhält man wieder aus der Kontinuitätsgleichung (1) 
ıst 


eu, 
27: F 1 
0 
also 
Die untere Integrationsgrenze ıst so gewählt, daß für y=0, v,=0ist. Für y: z RR 
", also 
4 
max A,x 7 d VE 
wonmt 
®, 
Un. 
1 max 4 


/ 
In Abb. 3 ist für !=0 und "=, aufgetragen. 
Max 
Die Ditferentialgleiehung für beliebige Werte von 7’ wird recht kompliziert, so daß ılır 
Lösung mit einem Rechenaufwand verbunden ist, der in keinem Verhältnis zum Ergebnis 
steht, da dureh die Interpolationsformel die Geschwindigkeitsverteilung für kleine Werte von 


!’ hinreichend genau angegeben ist. Die zweite Näherung für 1’ — pn gehoreht einer Di“ 


rentialgleiehung dritter Ordnung mit periodisch veränderlichen Koeffizienten. Im Falle 
wird die Differentialgleichung von solcher Form, daß sie nur auf graphischem Wege lös 
erscheint. Da in dieser Arbeit nur die erste Näherung für den Vergleich mit den Messung." 
in Betracht gekommen ist, so soll auf die Rechnung von weiteren Näherungen verziel 
werden. 


| Die zweite Näherung im Falle ? = ,,_ wird in einer Göttinger Dissertation von Herrn G. Cordes 


handelt werden. 


05 
RT SE | 
| A | | 
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3, Temperaturverteilung. a) Es seien sämtliche Stäbe erhitzt, sodaß die an- 
-smende Luft beim Durchgang durch das Gitter erwärmt wird. Die Wärmeausbreitung 
ster dem Gitter kann berechnet werden, sobald die Gesehwindigkeitsverteilung bekannt 
Die Gleiehung für den Wärmetransport ist‘) 


or Oy\ 


Da (ile 
D; »& den Temperaturunterschied irgendeines Punktes gegenüber der durehsehnittlichen Tem- 
ratur 9, bezeichnet, kinematische Austauschgröße. In Analogie zum Früheren kann 
ı Forms] für © ansetzen 


digkeit. Dim, 


‚o ©, die Temperatur des Punktes (4,0) angıbt. 
Damit wird 


oy\U dy 


ler integriert 


| 
N dd 


00 
\ Piıydy Uldy 
Yı 
war Die Austauschgröße z kann dadurch folgendermaßen ausgedrückt werden : 
U 
d d ), — 


Nach dem Prandtlschen Ansatz für die scheinbare Schubspannung kann Gl. (2) geschrieben 
werden: 


N) 
u ( ( (9). 
4 04 
wird IB... 
Durch denselben Ansatz wie früher für «,, den wir hier in der Form schreiben: 


ergibt sich auf genau dieselbe Weise wie oben 


Y 


( 
dy 


Setzen wir die beiden Größen x einander gleich, d. h. nehmen wir an, daß die Ausbreitung 

er Temperatur und der Geschwindigkeit durch den gleichen turbulenten Vermischungsvor- 

ang stattfindet, so erhält man folgende Beziehung zwischen der Verteilung von Geschwindig- 
laß ihr: seit und Temperatur 


roehnis 


Yı 
le —0 dy dy \dy/yzyı 
]Jösbar 

(y) 

ssung‘n Vählt man 9,—=0, so ergibt sich wegen =0, konst. 
Y, ‚so ergibt sich wegen F 


nach Prandtl ist die Temperatur- und Geschwindigkeitsverteilung bis auf Maßstabs- 


identisch. 
rdes 


Taylor: a.a.0. S. 01. 


ıls Intor. 29 
| 
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Nach Taylor lautet die Differentialgleiehung zur Berechnung der Störungsgeschwin: 
keit in erster Näherung 


N N? 
ou, 


wobei für 2 der Wert aus der Ableitung von Prandtl übernommen wird. Für u, setz 


wir an 


und erhalten damit wieder für die Austauschgröße 


dy 
Da 7° konst. ist, so unterscheiden sieh die Konstanten A, und A, nur durch den Faktor >. 
indem 1, 2.1, wie man dureh Einführen der Ansätze (10) und (12) in die Differentialgleiehungen 


bzw. (Ih) leicht verifiziert, die Geschwindigkeitsverteilungen nach Prandtl und Taylor 


werden identisch. Für "=, ist ,(y)= cos für wird f,(y) mit der bereits 


tabulierten Funktion F, identisch. Setzen wir wieder die Austauschgrößen einander gleich, 
so erhalten wır: 


Ist ft), d.h. die Gesehwindigkeitsverteilung bekannt, so kann man die Temperaturverteilung 

nach Taylor berechnen. Für wird die obige Gleichung durch = cos 2a 

/ 


erfüllt, also sind die Temperaturverteilungen für die beiden Annahmen in diesem Falle 
ilentisch. 


Für "=0 d.h. wird f” (0) d. h. die rechte Seite der Gleichung ver- 
schwindet, Da \ Pop dy nicht verschwindet, so muß (0) osein, d.h. P x, 


day, so gewählt werden kann, daß P’(y,) endlich ist. Eine Kurve, die für 9-0 eine senk- 
rechte Tangente hat, ist vom physikalischen Standpunkt aus natürlich unbrauchbar. Die 
lösung «der zugehörigen Differentialgleichung soll daher nicht weiter verfolgt werden. Diese 
senkrechte Tangente der Temperaturverteilungskurve rührt daher, daß in der Kurve der Ge- 
schwindigekeitsverteilung für der Krümmungsradius des Geschwindigkeits- 
profils an dieser Stelle also gleich Null wird. 

b) Der gzeheizte Einzelstab. Im Falle des gewöhnlichen Ansatzes für den Aus 
tausch 2 müßte die geheizte Zone hinter dem Einzelstab nicht breiter als „= + - nach 


jeder Seite werden, da der Austausch für 9 + „, wegen der horizontalen Tangente des G«- 


schwindigkeitsprofils verschwindet 0), Die Kurven der Temperaturverteilung sind dureh 
0% 


die Untersuchung unter a ermittelt. © bedeutet dann den Unterschied gegenüber der durel- 
schnittlichen Temperatur, ©, entsprechend für den Punkt (4,0). Die Messungen zeigen aber 
(siehe später im experimentellen Teil), daß eine wenn auch langsame Ausbreitung der Ten 
peratur nach der Seite stattfindet, je mehr man sieh vom Stab entfernt. 

Beim erweiterten Ansatz für den Austausch wird: 


(m) 


\ 


Ua 


also unabhängig von 7. Die Differentialgleichung für die Temperaturverteilung lautet \w 
vorher 


| 


— 


lo 


eits 


ich, 


lung 

/ 


itie 


enk- 

Die 
)jese 
‚eits- 


Aus- 


nach 


ureh 


ireh- 
ab 
Yen) 


\\ 
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»ezeiehnet b die halbe Breite des Temperaturnachlaufs, so muß, wenn von Wärmeverlust 
\uirch Strahlung abgesehen wird, die Beziehung gelten 
h 


\Ady=konst . 


\lan setzt nun an | 
9= fin) 

wo f(y) und zwei noch zu ermittelnde Funktionen sind und 7 = 4(r). Es ergibt sich 

dann für jeden Wert von 


h h 
dn = konst. 
h h 


/;ur Lösung der Differentialgleichung dient der obige Ansatz. Mit 
N 

0,04 


Fo) 


. 


wird die Differentialgleiehung 

die in folgende zwei Gleichungen gespalten werden kann 

gie. 

Die erste Gleichung gibt einmal integriert 

wo die Integrationskonstante gleich Null gesetzt wird, damit für 70, fo)>=0. Noch- 


malige Integration liefert 


InfmM)= — +In 
also 


Die zweite Gleichung liefert nach ‚Integration 


wo a eine Integrationskonstante ist, deren Bedeutung sich im folgenden ergibt. Damit wird®) 


‚\2 1 \ 
— f(V) exp I — 


Aa 


chreibt man ©, statt f(0), wo ©, die Temperatur an der Stelle (= — 2,11828.., he- 


sichnet, so wird 
3 / 
19 nd . 


In 


einem Schnitt == konst. ergibt sich somit die Gaußsche Wahrscheinlichkeitsfunktion 
1 die Temperaturverteilung. Für #=a ist 9-0 für jeden Wert von 9, d. h. die Wärme- 
Sbreitung geht von diesem Punkt aus, also ist die Achse des eeheizten Stabes im Punkte 
') anzubringen. Um die seitliche Ausbreitung der Wärme festzustellen, kann man nach den 


4 


y 
exp. statt „e hoch“ (z. B. exp. (— 2?) statt e =") 
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Punkten fragen, wo die Temperatur z. B. die Hälfte der Temperatur in der Mitte des gleich. 
Profils beträgt. Diese Punkte liegen auf der Kurve 


/ \1 
Y konst. (In 
Durch Änderung des Konstanten kann man eine Kurvenschar erhalten, wobei jede Kurve « 
konstantes Verhältnis zwischen einem Punkt der Kurve und dem Punkt auf der #-Achse |; 
gleichen Querschnitt konst. angibt. In Abb. ist die Ausbreitung der Temperatur angegehe 


| "25 
| | 
user | | | | 


Abb. 5. Temperaturverteilung in Schnitten x konst. für I!’ 


Die Temperaturprofile für «= = konst. ergeben sich dureh entsprechende Verzerrung der Gauß:- 
schen Wahrscheinlichkeitsfunktion, die in Abb. 5 für verschiedene Werte #: a aufgetragen ist '"). 


II. Experimenteller Teil. 


4. Versuchsanordnung. Die theoretischen Untersuchungen wurden durch Messungen im 
großen Windkanal der Aerodynamischen Versuchsanstalt in Göttingen nachgeprüft''). Die 
Messungen bezogen sieh dabei auf Geschwindigkeits- und Temperaturverteilung. Die Nach- 
laufströmung wurde mittels zylindrischer Stäbe aus Messingrohr von 1,5 em Durchmesser 


Abb. 6. Ansicht der Versuehseinrichtung. 


erzeugt. Die Stäbe waren 1,5 m lang, und oben und unten an waagerechten Balken befesti-' 
(vgl. Abb. 6), die wieder an einem großen Holzrahmen festgeschraubt wurden. Der Ho 
ralımen war am Fußboden und an der Bühne über der Düse des Windkanals mittels Schrau' 
zwingen befestigt. Die Messungen wurden in der Regel nachts ausgeführt, da die Versuc' 


10, Zahlenwerte nach E. Jahnke-F. Emde:a.a. O.S. 
11) Beschreibung des großen Windkanals siehe Ergebn. der Aerodyn. Versuchsanstalt zu Göttingen. I. Lief., > 
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‚nlage tagsüber für andere Arbeiten in Anspruch genommen war. Der ganze Aufbau mußte 
\aher nach jeder Schicht abmontiert werden. Der Rahmen mit len waagerechten Balken und 
nkreehten Stäben blieb während der ganzen Zeit, die die Messungen beanspruchten, in einem 
stück, so daß der Aufbau am Anfang und .der Abbau am Schluß jeder Sehieht nur geringe 
/eit in Anspruch nahm. Die Entfernung von der Ebene des Austrittquerschnitts der Düse 
is zur Ebene der Stäbe betrug 57 em, die Entfernung der Stäbe betrug zuerst 12, später 
; em (von Mitte zu Mitte gerechnet). Die Profile für Staudruck und Temperatur wurden in 
\orizontaler Riehtung in Höhe des waagerechten Durehmessers der Düsenöffnung gemessen, 

Abständen von 20 bis 125 em (hinter der Stabebene). Der Staudruck wurde mit Pitotrohr 
nd Manometer gemessen. Es wurde ein mit Alkohol gefülltes Vertikalmanometer Prandtl- 
‚cher Bauart benutzt. Das Pitotrohr wurde durch eine Lanze in den Luftstrom hineingeführt. 
Der Querschnitt der Lanze war hohl und von der Form eines aerodynamisch günstigen Streben- 
profils. Vom Pitotrohr führte eine Schlauchleitung durch den Hohlquerschnitt der Lanze 
sum Manometer. Die Lanze konnte mit Hilfe einer Verschiebeeinriehtung in horizontaler und 
vertikaler Riehtung bewegt werden, so daß jeder Punkt einer Meßebene auf diese Weise er- 
reicht wurde. Für kissenazen in verschiedenen Ebenen mußte die ganze Verschiebeeinrichtung 
in Riehtung der Strahlachse verschoben und justiert werden, damit sich die Öffnung der 
Sonden in der gewünschten Entfernung vom Gitter parallel zur Düsenkante bewegte. 
Der statische Druck wurde mit Hilfe einer statischen Sonde gemessen. Da sich heraus- 
stellte, daß dieser weniger als 1°, des Gesamtdruckes betrug, und da aus früheren 
theoretischen Untersuchungen hervorgeht, daß er in der Tat noch kleiner ist, so wurde er 
vernachlässigt'?). Wie zu erwarten war, wies der gemessene statische Druck Schwankungen 


auf, die für y= + , “und 4=0 besonders groß waren (n=1,3,5 usw.). Die Windge- 


schwindigkeit betrug rd. 20 m/sek. 
Die Temperaturmessung erfolgte durch ein Thermoelement aus Konstantan und Kupfer. 
Die eine Lötstelle wurde in den Luftstrom hinter den geheizten Stäben, die andere Lötstelle 
den ungeheizten Luftstrom gebracht (Abb. 7). Die thermoelektrische Spannung wurde 


I =lötstelle 7 
V = Voltmeter 
o=geheizter Stab (5.9) 
V o=ungeheizter ” (1.410..3) 


Abb. 7. Anordnung zur Messung der Temperatur (schematisch). 


mittels eines empfindlichen Zeigergalvanometers mit Fadenaufhängung von Siemens u. Halske 
gemessen. Auf dem Spiegel ist das Bild des Zeigers angebracht, das auf die Instrumentskala 
projiziert wird. Die Trägheit ist infolge dieser Anordnung sehr klein und die Einstellzeit 
sehr kurz (3 Se kunden). Das Thermoelement wurde im Anschluß an die Messungen durch 
Kintauchen der einen Lötstelle in dünnem Öl geeicht. Die Temperatur des Öls wurde in 
einem Dewar-Gefäß konstant gehalten. Die Öltemperatur wurde durch Hinzugießen von 
heißem Ol geändert und die neue Temperatur mit Hilfe eines Normalqueeksilberthermometers 
bestimmt. Die Temperatur der anderen Lötstelle wurde in einem Wasserbad nahezu kon- 
stant gehalten 

Nur die fünf mittleren Stäbe (Stab 5 bis 9) wurden elektrisch geheizt, da diese Anzahl 
für die Ausbildung eines periodischen Verlaufs der Temperatur genügte. Diese Stäbe waren 
wie die anderen hohle Messingrohre von 1,5 em Durchmesser, durch die Rohre von 0,8 em 
Durchmesser hindurchgesteekt waren. Um das innere Rohr waren Drähte gewickelt, die 
durch aufgezogene Glasperlen gegen die Rohrwände isoliert waren. Durch den Draht eines 
liohres wurde. ein Strom von 2 Ampere geschickt. 


3. Versuchsdurchführung Zuerst wurden die Druckmessungen durchgeführt. Der Ab- 
‘sand der Stäbe betrug dabei 12 em, später 6 em. Es wurde in 30, 50, 75, 100 und 125 em 
"ntfernung gemessen, die Druckprofile sind in Abb. S aufgetragen. Die Auswertung der 
\lessungen ergab, daß die Gesch w indigkeitsprofile erst von 75 em an so ausgebildet waren, 
vie es die Theorie voraussetzt. Für die Entfer nung e=30 und e=50 em war das Geschwin- 
ekeitsprofil in der Mitte zwischen zwei Stäben zu flach, gerade hinter den Stäben zu spitz. 

\ährend jeder Arbeitsschieht wurde die Temperatur der Luft und des Alkohols, sowie der 
"rometerstand abgelesen. Temperaturmessungen wurden nur für einen Abstand der Stäbe 


=), W. Tollmien: ZAMM, Bd. 6 (1926), S. 468. 
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gleich 6 em ausgeführt, weil bei =12 em in der Entfernung, die einem sinusförmig au 
zebildeten Profil entsprach, keine meßbaren Unterschiede in der Temperatur bestanden. | 
zeigte sich nämlich, daß die Wärme so schnell ausgetauscht wurde, daß im Gebiet hinter d: 
geheizten Stäben bereits 50 em von der Stabebene entfernt keine meßbaren Temperaturunt: 
schiede bestanden. Dagegen konnte der Temperaturverlauf hinter dem geheizten Einzelst: 
in «> em Entfernung noch gut verfolgt werden. Im Gegensatz zur Geschwindigkeitsmessun: 
die wegen der Schwankungen im Manometer nach Einstellung eines neuen Meßpunktes zien 
lich lange Zeit in Anspruch nahm, gingen die Temperaturmessungen rascher vonstatten. (; 
das Voltmeter sich sehr schnell einstellte. 
6. Auswertung des Versuchsmaterials. 


für A=12 em sind in 


Die Gesamtdruckprofile 


Abb. S wiedergegeben. Da der Druck in der Düsenvorkammer mit Hilfe der Druckwaaso 
P 
| 
N 
2 N 4 | 


0 
A 


hinter dem Gitter für 
erößte Druck in der be- 
treffenden Ebene), 


Abb. s, Gesamtdruekprofile 
/ !em ist der 


Abb. 9. Gesamtdruekprofile hinter dem Gitter für 


A=6e 
bem. 


während einer Meßreihe konstant gehalten werden konnte'”), so wurde das Meßprotokoll naeh 
nebenstehendem Schema ausgewertet, wo als Beispiel e= 125 em gewählt ist und die Drücke 
p in Millimeter Alkoholsäule angegeben sind. Außerdem sind die Gesamtdruckprofile für 


tee) 
Tabelle 2. 


p 
Hnax 
SS 29,75 0,974 0,087 
30,35 0,994 0.997 
05 30,55 1.000 1.000 


6b em und kleinere Werte e in 


wurde, 
mußte. In 
und 


Abb 


em aufeetraeen. 


Abb. 9 dargestellt. 


13) Siehe die unter 11 angegebene Quelle, S. 20. 


Die Festlegung der Koordinatenachlıs: 
4=0 geschah in der Weise, daß die Symmetrielinie im Profil hinter jedem Stab gezoge 


Als Kontrolle ergab sieh, daß die Entfernung dieser Symmetrielinien gleich 4 sei 
und 11 sind die entsprechenden Geschwindigkeitsproflle für 
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N Abb. 10. Geschwindigkeitsprofile hinter dem Gitter Abb. 11. Gesehwindigkeitsprofile hinter dem Gitter R 
für 12 em. für 2 


Die Meßprotokolle für die Temperaturverteilung wurden in entsprechender Weise wie ' 
N bei der Druckmessung ausgewertet. Die Ergebnisse der Messungen sind für den Fall, daß 
> Stäbe geheizt wurden, in Abb. 12 bis zur Mitte von Stab 7 (Symmetrieachse) aufgetragen. 
Die Temperatur wird demnach erheblich schneller ausgetauscht als die Geschwindigkeit, da 
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» 12. Temperaturprofile hinter vielen geheizten Abb. 13. Temperaturprofile hinter dem 
Stäben. geheizten Einzelstab (b’ bezeichnet die Breite 


in halber Höhe). 
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bereits in 50 em Entfernung hinter dem Gitter die ursprünglich wellenförmige Temperatı 
verteilung in eine „Hochebene“ ausgeartet ist. Für den geheizten Einzelstab sind die Me 
ergebnisse in Abb. 13 aufgetragen. In diesem Falle ließ sich die Temperaturverteilung h 
sonders gut verfolgen, so daß man in größerer Entfernung als 50 em noch gut verwertbaı 
Messungen erhalten konnte. Die Eiehung des Thermoelements ergab, daß ein Skalenteil a 
Amperemeter 1/7° (genauer 0,142°) Celsius entspricht. Das Ergebnis der Eichung ist in Tabelle 
wiedergegeben. Dabei bedeuten #, und 't, die Teinperatur der kalten bzw. warmen Lötstell 
s die angezeigten Skalenteile des Voltmeters. 


Tabelle. 


16,4 24.6 8.2 0,142 
16.4 24,2 55,1 7.8 0,142 
16.4 23.4 49,3 7.0 0.142 
16.45 22.6 43,0 6,15 0,143 
16,5 22,6 12,8 6.1 0.143 
16.7 22,1 37,2 0,145 
16.8 21,35 32.5 4.55 0.140 
16.9 20,2 23.4 3,3 0,141 
16.0 19,4 17.8 2,5 0.141 
16.0 IN.6 12.0 1.7 0.142 


Im Mittel 42:5 = 0.142 


7. Vergleich der Geschwindigkeitsmessung mit der Theorie. a) Das Potenzgesetz. 
Es soll zunächst nachgeprüft werden, ob die Abnahme der Maximalgeschwindigkeit mit der 
ersten Potenz der Entfernung «= richtig ist. Die Formel für die Störungsgeschwindigkeit u 


lautet für 7 (9) 


In 
- 
Sa’\l /. 
Für ist u, = U... — U, also 
konst 


wobei in der Konstante der Mischungsweg I! als Unbekannte enthalten ist. Wertet man «die 
Auftragung der Abb. I1 entsprechend aus, so ergibt sich die Tabelle 4. 


Tabelle 4. 


e (em) max 
50 0,025 
0.943 0.0606 
100 0,957 0,0449 
125 0,065 0,0363 


Man trägt aus Tabelle 4 die Werte Tr als Funktion von x auf, wobei x mit 


nicht identisch zu sein braucht, da der Nullpunkt des Koordinatensystems noch nicht fe 
gelegt ist. Für die drei größeren Entfernungen (75, 100 und 125 em) ergibt die Auftragu' 
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\bb. 14), daß die Punkte ziemlich genau auf einer Geraden mit der Neigung 1:1 liegen, 
I. das Potenzgesetz ist erfüllt und . ist gleich e zu setzen, da eine Verschiebung der 
‚nkte um gleich große Strecken nieht notwendig ist. Der Vergleich der Temperaturmessung 
‘ter dem Einzelstab mit der Theorie läßt auch darauf schließen, daß der Nullpunkt des 

‘sordinatensystems in kurzer Entfernung vor dem Gitter liegt. Für e=50 em ist das Po- 
nzeesetz nieht mehr erfüllt, d. h. die erste Näherung trifft für diese Entfernung nicht mehr 

In Abb. 15 ist die Auftragung in gewöhnlichem Maßstab gemacht, damit man besser den 

‘dl der Genauigkeit sehen kann, mit welcher das Potenzgesetz erfüllt ist. Der Geraden 

‚ıf Logarithmenpapier entspricht dabei eine gleichseitige Hyperbel in gewöhnlicher Auf- 


(ragung. 
00: 
>| 005 
> 
403 
> 
0.02 SID 
S | 
0,07 0 | | | 
20 cm 20 25 75 700 725 cm 
— 
Abb. 4. Nachprüfung des Potenzgesetzes in log- Abb. 15. Nachprüfung des Potenz- 
arithmischer Auftragung. gesetzes in gewöhnlicher Auftragung. 


b) Die Mischungswegkonstante j: Nach Festlegung des Koordinatenanfangs- 
punktes sind wir nunmehr in der Lage, die einzige empirische Konstante der "Theorie, näm- 
lieh die Misehungswegkonstante zu ermitteln. 


Nach Gl. (3) ist 


U [A\®A 


also für y=0 
U | 


Da in dieser Formel nur T unbekannt ist, so kann ermittelt werden. Das Ergebnis ist in 


labelle 5 zusammengestellt. 


Tabelle». 


r (em) 8a 
75 0.0606 1.550: 109 0,103 
100 0.0449 2,067 103 0,104 
125 0.0363 2,583. 109 0.103 
Ist also im Mittel 1=0,1034 für 7’ Im Falle ist nach Früherem Gl. (7) 


1,161 0,1204. 
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Nimmt man statt / als charakteristische Länge die Breite des Nachlaufprofils in halber Hö 
so ist diese im vorliegenden Fall gleieh der halben Gitterentfernung zu setzen, also nach « 


von Schliechting'*) benutzten Bezeichnungsweise 25’ = ,. Damit wird I,—=0,48b’; währe: 


nach den Schliehtingsehen Versuchen 1,-=0,47 5b’ herauskommt. Die Nachlaufmessung 
hinter dem Einzelstab werden dureh diese Übereinstimmung in ihrer Richtigkeit bestäti: 


e) Vergleich der Messungen mit der Rechnung in Schnitten z=konst. 
man in Abb. 2 die Meßpunkte für die fünf Entfernungen 30, 50, 75, 100 und 125 em ein, so ergi)i 
sich Abb. 16. Man sieht, daß die Meßpunkte für die drei größeren Entfernungen in guter 


e=-T25cm 
oe=-75 " 
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Abb. 16. Vergleich zwischen Theorie und Versuch bei der Geschwindigkeitsverteilung. 


Übereinstimmung mit der Theorie sind, während für e = 30 und e=- 50 em die Meßpunkte ziem- 
lich herausfallen. Dies liegt daran, daß die Theorie der ersten Näherung für diese kleineren Ent- 
fernungen nieht mehr zutrifft, wie es auch mit dem Potenzgesetz der Fall war (vgl. unter a). 
Aus Abb. 16 ıst ferner zu ersehen, daß die Streuung der Versuchspunkte auch für größere 
Entfernungen zu stark ist. um auf einen genaueren Wert für U’ schließen zu lassen. Da je 


doch die Versuchspunkte zwischen den beiden für U =-0 und I’ Ri geltenden Kurven liegen, so 
ist, wie man bisher auch annahm, mit einem kleinen Wert für !’ zu rechnen. Die Annahme 
die der Aufstellung der Interpolationsformel (5) zugrunde gelegt wurde, erscheint 
somit berechtigt. 


8. Vergleich der Temperaturmessung mit der Theorie. Die Temperaturmessungen wurden 
für 16 em ausgeführt. Da die Temperaturverteilung nach dem Ansatz von Prandtl mit 
der Geschwindigkeitsverteilung identiseh sein soll, während nach der Annahme von Taylor 
(die Geschwindigkeitsverteilung bekannt sein muß, um die Temperaturverteilung zu berechnen. 
war es erforderlich, die entsprechenden Geschwindigkeitsmessungen auch für 1=6 em aus 
zuführen. Die Ergebnisse dieser Messungen sind bereits in Abb. 9 und 11 wiedergegeben. 


Kin Vergleich mit der Theorie für viele geheizte Stäbe ließ sich nicht in dem Malw 
durehführen, wie bei der Gesehwindigkeitsverteilung. Wie bereits früher erwähnt, liegt «dıvs 
daran, daß die Wärme viel schneller als die Geschwindigkeit ausgetauscht wird. In kurzer 
Entfernung hinter dem Gitter, wo gut meßbare Temperaturunterschiede vorkommen, ist dus 
Temperaturprofil in Schnitten «= = konst. mieht genügend ausgebildet. In größerer Entfernunz. 
wo die Annahmen der Theorie gut erfüllt sind, kommen kaum noch meßbare Temperatu' 
differenzen vor. Um eine Vorstellung von der Temperaturverteilung in großem Maßstab 
geben, sind in Abb. 17 die theoretischen Kurven mit den Meßpunkten für e gleich 15, 
20 und 25 em eingetragen. Die Meßpunkte schließen sich nur in großen Zügen an die tlıc 
retischen Kurven an. 
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Für einen geheizten Stab im Gitter (Stab 7) ergab die theoretische Lösung mit 


die Gaußsche Wahrscheinliehkeitsfunktion für die Temperaturverteilung in Schnitten 
konst. Die Kurven der Abb. 13 lassen auch einen solehen Verlauf vermuten. Trägt man 


‘» Versuchspunkte für die Entfernungen e= 30, 37,5, 45 und 50 em so auf, daß die Kurven 


Abb. 17. Vergleich zwischen Theorie und Versueh bei der Temperaturverteilung. 


y=0, 9= 9, und zusammenfallen, so ergibt sich in Abb. 18, daß 
die Versuchspunkte auf einer ziemlich engen Straße liegen. Normiert man die theoretische 
Kurve auch in der Weise, daß sie in und y init den Versuchspunkten zu- 


sammenfällt, so wird die Übereinstimmung zwischen Theorie und Versuch befriedigend mit 
Ausnahme der Punkte, die am Rand der erwärmten Zone liegen. 


Anschließend wurde die Temperaturverteilung hinter dem geheizten Einzelstab ohne 
Nachbarstäbe gemessen, wobei die Messungen von Fage und Falkner als richtig bestätigt 
werden, die zur Prüfung der Rechnung von Taylor für den Einzelstab ausgeführt wurden '®). 
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'S,. Vergleieh zwisehen Theorie und Versuch bei der Temperatur- Abb. 19. Vergleich zwischen Theorie und 
lung hinter dem Einzelstab, e ist diejenige Kurve, die sieh den Versuchbeider Temperaturverteilung hinter 
ichspunkten am engsten anschließt, b die theoretische Kurve dem geheizten Kinzelstab ohne Nachbarstäbe 
(Gaußsche Wahrscheinlichkeitsfunktion). ist die ıiheoretische Kurve nach G, 


l. Taylor, b und e haben dieselbe Be- 
deutung wie in Abb, 18). 


15), A. Fage und V. M. Falkner: Proe. Roy. Soe. A., Vol. 135 (1932), S. 702, 
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Auch diese Punkte sind in Abb. 15 eingezeichnet. In Abb. 19 sind nochmals die gemesseı 
Punkte eingetragen, ferner die theoretische Kurve nach Tavlor (a) sowie die Gaußsc 
Wahrscheinlichkeitsfunktion (5b). Gestrichelt eingetragen ist diejenige Kurve, die sich 
gemessenen Punkten beim geheizten Einzelstab im Gitter am besten anschließt (ce). Dal 
sind die Kurven in derselben Weise wie früher normiert. Die wirklichen Breiten des Na 
laufprofils sind für den Einzelstab größer ohne als mit Nachbarstäben. Die in Abb. 19 wie 
gegebene Messung bezieht sich auf die Entfernung e=30 em, wobei eine Breite des Na: 


laufprofils in halber Höhe (für Oo = 


> |gleich 63 mm, während beim Gitter entsprechend ı 


Breite 55 mm gemessen wurde. Dieser Unterschied rührt daher, daß beim Gitter der A 


tausch für + wegen sehr klein wird, während beim Einzelstab ohne Nacı 


barstäbe noch yr -=0. Es mag bemerkt werden, daß bei der Berechnung der Temperatur 


verteilung hinter einem rotationssymmetrischen Körper eine nur von .« abhängige Annalıme 
über den Austausch für die Temperaturverteilungskurve ebenfalls die Gaußsche Wahrsehein- 
licehkeitsfunktion liefert. Die Übereinstimmung mit der Messung ist sehr gut, wie bisher un- 
veröffentlichte Göttinger Versuche zeigen. Allerdings fehlt im rotationssymmetrischen Fall 
für diese Annahme über den Austausch die theoretische Begründung. 


Zusammenfassung. Mit dem Ansatz von Prandtl für die scheinbare Schubspannung 
wird die turbulente Geschwindigkeitsverteilung hinter einem Gitter in erster Näherung be 
stimmt. Dabei wird von dem erweiterten Prandtlschen Ansatz für den Austausch Gebrauch 
gemacht. In diesem Fall ist die Geschwindigkeit sinusförmig verteilt. Mit dem gewöhnlichen 
Ansatz für den Austausch ergibt sich eine davon etwas abweichende Geschwindigekeitsver: 
teilung. Die maximale Nachlaufgeschwindigkeit nimmt proportional der Entfernung vom 
Gitter ab. Weiter wird die Temperaturverteilung nach dem Ansatz von Taylor unter Be 
nutzung des Prandtlschen Ansatzes für den Austausch berechnet, wobei für den erweiterten 
Ansatz «die Temperaturverteilung ebenfalls sinusförmig wird. Mit dem vereinfachten Ansatz 
bekommt die Kurve der Temperaturverteilung nach Taylor an den Stellen der Temperatur- 
maxima und -minima senkrechte Tangenten. Diese eigentümliche Singularität hängt mit dem 
vereinfachten Ansatz für den Austausch zusammen. Für einen geheizten Stab im Gitter 
ergibt sich im Falle "= R- die Gaußsche Wahrscheinlichkeitsfunktion für die Temperatur- 


verteilune. 


Die Versuche über die Geschwindigekeitsverteilung bestätigen die bei der Rechnung der 
ersten Näherung gemachte Annahme für Entfernungen etwa gleich dem sechsfachen «es 
Gitterabstandes. Für kleinere Entfernungen ist das Geschwindigkeitsprofil hinter den Stäben 
zu spitz, in der Mitte zwischen zwei Stäben zu Nlach. Die Messungen der 'Temperaturver- 
teilung hinter vielen geheizten Stäben waren wegen des schnell erfolgenden Austausches der 
Temperatur zu einem Vergleich mit der Theorie wenig geeignet, dagegen war die gemessene 
Temperaturverteilung hinter einem geheizten Stab mit ungeheizten Nachbarstäben in grobe 
Entfernung hinter dem Gitter in guter Übereinstimmung mit der Theorie. Die Temperatur 
verteilung hinter dem zeheizten Einzelstab wurde ebenfalls gemessen und in Übereinstimmung 
mit früheren englischen Messungen gefunden. Die Mischunesweerkonstante 1 ergibt sieh 
fast genau so groß wie bei der Nachlaufmessung hinter dem Einzelstab 0,124, 4 = Gitt‘ 
entfernung). 


Die Anreeune zu dieser Arbeit verdanke ich Herrn Professor Dr. L. Prandtl. 
mir auch Gelegenheit gab, die Versuche in der Aerodynamischen Versuchsanstalt zu Göttingen 
auszuführen. Herr Dr. H. Scehliehting (jetzt in Friedrichshafen) hat während meiner .\" 
wesenheit von Göttingen für die Fertigstellung der Versuchseinrichtung Sorge getragen un 
auch während der Versuche diesen sein Interesse gezeigt. Herr W. Frössel und Herr 
Hedewig haben mir während der Versuchsdurchführung wertvolle Hilfe geleistet. Den 
nannten Herren sei auch an dieser Stelle bestens gedankt. Für finanzielle Unterstütze 
bin ich Norsk Varekrigsforsikrings Fond zu Dank verpflichtet. OL 
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hene und räumliche Strömung zäher, inkompressibler, 
ägheitsfreier Flüssigkeiten zwischen exzentrischen, relativ 
zueinander rotierenden Zylinderflächen. 
Beitrag zur Theorie der Schmierung. 
Von H. Reissner in Berlin-Charlottenburg. 
Mitteilung 
B. Räumliche, zähe Strömung im Ringraum. 


1. Vorbemerkung. In Mitt. I wurde das ebene Sommerfeldsche Problem der zähen 


-trömung im engen Ringspalt exakt gelöst. Für das räumliche Problem lag bisher nur die 


\rbeit von A. G. M. Michell und die an sie anschließenden Arbeiten über die zerade 
\eilströmung vor. In der folgenden Mitteilung gelingt es auch die räumliche Aufgabe für 
Kinzströmung zu lösen. 


Wie in Mitteilung I werde eine Bezugslänge und zwar der Radius a des inneren 
ireises (Abb. 1) eine Bezugsgeschwindigkeit 4 und zwar die Geschwindigkeit am inneren 
Kreise und ein Bezugsdruck 7" (wo 7 die Zähigkeitskonstante) eingeführt, so daß als dimen- 

Aa / 5 


sionslose Variable wie in Mitteilung | auftreten: 


Die Gleichungen der stationären, inkompressiblen Stimmung lauten dann wiederum: 


2. Integration durch Partikulärintegrale. Nach dem Vorgange von A. G.M. Michell bei 
der Lösung für den von Ebenen begrenzten Schmierkeil befriedige man zunächst die Be- 
dinzung, daß die Normalspannungen 


| Om 
22) ) 
UT 
an den Enden z >». Null sein sollen. 
a 2a 
Man setze mit 
 Pn C0Sa„?, N -w.sma„2 . 


Die Winkelfunktionen für die Geschwindigkeiten sind hier so angesetzt, wie es der 
\ulbau der Gl. (2) erfordert. 


Die Gln. (2) nehmen dann bei Einführung des Operators 


2 
Ay 


folgende Form an?): 


p Ph 


Tau Wu = sc). 
') Mitteilung I in Bd. 15 (1955), H. I m. 2, 8. 81 bis 87 dieser Zeitschrift. 
°) Dieses Gleiehungssystem kann man auch so schreiben, daß in jeder Gleichung nur eine Variable auftritt, nämlich 


ı durch Besselsche Funktionen integrieren lassen. Jedoch wird von dieser Form im folgenden kein Gebraueh 
ht, 


15* 


Sse] 
h d 
Da) 
vier 
nd die 
r 
| 
| | 
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3. Einführung von Polarkoordinaten. Nach Abb. 1 führe man Polarkoordinaten « 
so daß für die Geschwindigkeitskomponenten zu setzen ist: 


Dann nimmt das Gleichungssystem (2e) die folgende Form an: 


r or 00° 
Io l Oug 

_2o 

Zofz- 


Abb. 1. 


In folgendem werde der volle Ringraum betrachtet und wegen der Einwertigkeit aller 
Geschwindigkeiten und Drucke und entsprechend dem Bau der GIn. (2d) gesetzt°): 


w„=2W,„xrsin(k d) 2, =ZP,xsin(k 9) 
Führt man noch Geschwindigkeitsaggregate durch die Ansätze 


Ul,n,k Us, n,k Us, n,k 
ka) 


ein, so nehmen die Gln. (2d) nach Einführung des Besselschen Operators 


Bi (a + 


or 


und bei vorübergehender Weglassung des Index n die folgende Gestalt an: 


d P, r d ıD 
r dr r r dr 


Zu beachten ist, daß die letzte der Gln. (2e) eine einschränkende Bedingung für 
Konstanten der Lösungen der ersten drei Gleichungen vorstellt. 


3 


Dieser Ansatz ist mit nicht notwendig ganzem k auch brauchbar, wenn an den Enden eines Bogens Yo : 
bei nieht ausgefülltem Ringraum die Spannungen an den freien Oberflächen gegeben werden. 


| 1 0?’ 


(2d), 


t aller 


(ta) 
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4. Die Koeffizienten der trigonometrischen Doppelreihen in (z, )) als Besselsche Funktionen 
s Radius r. Man kann die Lösung nun in folgenden Schritten gewinnen: 


Für #0 folgt aus den Gln. (2d) und (4), daß nur 7’, „,o auftritt, und zwar ergibt sich 
ior wie in folgendem bedeutet wie üblich 
AJ,(tar)+ BN,(iar) 


‚ mit zwei Integrationskonstanten behaftete Aggregat der Besselschen Funktion imaginären 
‚reuments k-ter Ordnung. 


k+0. 
ür ergeben die Gln. (2e) 
B,(W,)= taZ,ltar) 
und wenn man die Differentialbeziehung 
beachtet 


ie weitere Integration wurde zuerst nach der Methode der Variation der Konstanten durch- 
veführt, lieferte aber auf diese Weise auch unter Ausnutzung aller Besonderheiten der Bessel- 
schen Funktionen sehr unhandlicehe Ausdrücke. 


5. Entwicklung der Lösung vom Punkte 7—1 aus. Da das betrachtete Strömungsgebiet 
nur von r—1lbis r=1-+0(0o 1) geht, ist es zweekmäßig, die auftretenden Besselschen 
Funktionen und auch die infolge der Störungsfunktionen der Gln. (2f) entstehenden partikulären 
Integrale als Potenzreihen der unabhängigen Variabeln o darzustellen, wobei man wegen der 
leinheit von o nur bis zur zweiten, höchstens dritten Potenz zu gehen braucht. 

Zunächst seien die Potenzentwieklungen für die homogen gemachten (von den Störungs- 
funktionen der rechten Seiten befreiten) GIn. (2f) angeschrieben, die man bequemer nach der 
\ethode der unbestimmten Koeffizienten, als aus der Reihenentwieklung der Besselschen 
Funktionen erhält. Diese Anteile der Lösungen seien mit dem oberen Index (0) bezeichnet. 
Die homogenen Lösungsanteile sind: 


Fir 
ur = 


» (u) a 4 | 
Für 7==0 
a@+(k — 1? a? +-3(k—1)% | 
0 a®+k— 1’ +2 
> (ba). 
2 ( k 1 + > 6; 0 + .. 
+2 
(0 (0) 
1? 4-2 


[3 

| 
We). 
ür 
Io 
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Die Partikulärlösungen für die Störungsfunktionen ergeben 
vereleiehung auf Grund der oben angeschriebenen Reihen von P,, 
diese Anteile bezeiehnen möge, wie folgt: 


> 


=* (k + 


6 


(0) 0 
1a 4, 


a A, 


Wie sehon oben begründet, 


sich dureh Koeffizient: 


wobei der obere Index 


A) 


-2) A, ”) (6 


enthalten diese Lösungen noch zu viele Konstanten, da =. 
statt der Kontinuitätseleiehung dive=0 die allgemeinere Bedingune I div ec 


0 erfüllen. 


Man muß also von den gegebenen Lösungen (6a) und (6b) noch verlangen, daß si 


identisch die Kontinuitätsgleichung, die letzte der 


>e) erfüllen. 


Benchtet man nun, daß, wie im nächsten Abschnitt sofort zu sehen ist, die innere 
Randbedingzung verlangt, daß 
werden, so führt die Kontinuitätseleiehung zu den folgenden Bedinzeungen: 


6. Die Befriedigung der Randbedingungen. 
und (6a) (6b) für 
der Kontinuitätsgleiehung 


bei 1% 


für jedes n 


0 zwei Konstanten und zwar und 


Das System der Partikularintegrale (5a) für 
-=0 enthält nach Beachtung der Koeffizienteneinschränkungen (7) 


bei fünf Konstanten, da von den acht Konstanten A®, BV, 09,09, 


drei wegen der Einsehränkungsbedingungen (7) und (9) eliminiert werden können 


Diese Konstanten müssen nun entsprechend den Randbedingungen für u,, ug und 


am Innen- und Außenrande bestimmt werden. 


Es ist folgendes zu verlangen: 


6a. Am Innenrande 7 = 1(o 0) 
uy = 1. 
Dies liefert für 
und da am Innenrande die Umfangsgesehwindigkeit | 
für die die Bestimmungsgleichung: 
\ 2 cv) | 
cos a,2 l |z von a; bis + 


und damit: 


= —(—1}' 


Ferner ergibt sich bei = 0 und k=E0: 


Aus ( Usn,k), „= 0 x” —(), 


womit die Voraussage von Gl. (9) bestätigt ist. 


6b. Am Außenrande (o 
am Außenrande (0 
und für 


nur eime Konstante 


sein sollte, so ergibt sieh 


Es bleiben demnach für die Erfüllung der Bedingung‘ 
o(#)) für nur noch drei Konstanten etwa und 
verfügbar und zu bestimmen, und man muß folgen‘ 


Funktionen als Geschwindigkeitskomponenten am Außenrande betrachten?). 


ı) Für Au,aP ist nach Gl. (7) 


gesetzt. 


mu 
| 


(6 


Nnere 


Für 
n (7) 


nnen 


nd 


t sich 


(). 


gen! 
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Je 


- 


U B 1,0 


6 
k +2 


> 6 


+ (Bu,k'” k — a) + 6 


+2. 
4 
® 
gar...) 


Die Randbedingungen fordern nun, daß am Außenrande (o= o,)",, ty und w ver- 
schwinden, daß also 


| 0 


‚snk9 (u). 


Beachtet man, daß die verhältnismäßige Lichtweite 0,(9) sehr klein bei dem hier vor- 
iezenden Problem, nämlieh <O,OL, sein soll, so erhält man aus den vorhergehenden Gln. (11) 
unter Berücksichtigung der Beziehungen (4a) zwischen T/,, U, und U,, U’, und wenn man 
noch zur Abkürzung = Br, Aus = Ar, = Dr setzt: 


\ Ssınk [B,k— D,al=0, \ cos Ar 0, 
> 
(lila), 
2 


wo nach Gl. (8) B,"” = Bu, : (— die Erregung der Strömung dureh die 


1)" 4 

(2n+1)a 
nlangsgeschwindigkeit = am inneren Rande wiedergibt und, wie es sein muß, mit 
3 V alle Konstanten ebenfalls Null werden würden. 


Da aus der ersten der obigen Gleichungen gefolgert werden kann, daß 


0 lassen sich die GlIn. (11a) auch vereinfachen zu: 


( Ya | n | 


\ sın (k 9) k A}; 


Mi 
B 
(7) 
a 
+...) 
| 
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In diesen Gleichungen ist zu beachten, daß o„ eine Funktion des Winkels #9 ist, u 
zwar bei sehr kleiner Exzentrizität e und sehr kleiner Radiendifferenz d (mittlerer Liehtwei: 
eesetzt werden kann: 


wo 
E 
oder mit 5 | 


wobei aus Symmetriegründen vorausgesetzt werden kann und übrigens auch aus der inneren 
Randbedineunge #, 0 folet, daß die größte bzw. kleinste Liehtweite an der Stelle 9- 0 
bzw. 9 — a liegt. 

Die Gln. (11b) schreiben sich nun bequemer, wenn man «die Konstanten A, und 7, 
dureh die folgenden Ausdrücke ersetzt: 


0 


Sie nehmen dann die folgende etwas übersichtlichere Gestalt an: 
4 d) (1 +yeosd) Isinkd|b- (1le), 
oder wenn man die linken Seiten in trigonometrische Reihen umformt: 


k+1 


| 
b, + ax 5 2 A 


(id). 


7 
sınk IR a’\dk— TAr +17) j 
Kliminiert man 5, mit Hilfe der zweiten dieser Gleichungen und führt man noch als ab- 
kürzende Bezeichnung ein 


(14), 


dann kann man als Rekursionsformel oder Differenzengleichung schreiben: 


cos lan ] I}; 1/7 + a7, (Ir ı + -] -4- 7 «os (lle). 


7. Die Lösung der Rekursionsformeln. Entwicklung nach Potenzen der verhältnismäßigen 
Exzentrizität „. Man konnte zunächst daran denken, auch die rechte Seite dieser Gleichung 
in eine trigonometrische Reihe zu entwiekeln und die entstehende inkomogene Differenzen- 
eleichung mit veränderlichen Koeffizienten zu integrieren’). 


Dieser grundsätzlich durchaus gangbare Weg erwies sich praktisch numerisch als unzweck- 
mäßig, jedoch führte die Entwieklung der gesuchten Koeffizienten a, nach Potenzen des 
Parameters z(y<1), d.h. nach der verhältnismäßigen Exzentrizität zu schrittweise und ein- 
fach lösbaren Bestimmungsgeleiehungen, wie folet: 


ls werde der Ansatz gemacht: 


Dieser Ansatz führte zu folgendem: Es werde gesetzt bzw. bereehnet: 


\ ( I)k 
T | /vos.f \ y2 7 + / \ j 


dann geht die letzte Gl. (1le) über in die Differenzengleichung: 


+ 1}. 
k+l 


Man kann für große A eine asymptotische Lösung mit an angeben und a, daraus bestimmen, 


hält jedoch sehr unbeqnueme Formeln. 
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2. »I dann die rechte Seite der Gl. (lle) nach dem binomischen Satz entwickelt, nach Potenzen 
wei ı y geordnet und der Faktor jeder Potenz von y für sieh Null gesetzt, so erhält man: 
( 


»se Gleichungen lassen sich nun schrittweise wie folgt lösen, wobei sich herausstellt, daß 


lesmal 
(19). 
für kpi, 
neren 
as eben die Berechnung besonders erleichtert ®). 
für alle 
(13) 
le) 0 N 
2k 
coskd [ar — + (Ir ı + ara 1,1 la] = (1 cos2n), 
(ld). 0 2 I, 2 
für k>2 
Is ab- 
| \' | 
9 3 | 3, 21i,(,+21)\ 
> | I, 2(,+21) 
chung 
ANZEn- für 
weck- \' 8 1 | 
n «des cosk far — 1,3 -ı + ar, -ı + + ar + 1,3 +5 cos 4 
( 1, A +1, +1, +1, 
2 +1) +41, (, +21) 


+1) (1, +1,) 


a,,=d, 


ars für k>4 usw. 


Ad 


. . . . 
‚ Im folgenden ist, wie aus einer bekannten, allgemeinen Formel sieh ergibt, beachtet, daß eos? 4 (1+00s29#), 


COS 4 cost cos29+ cos4%$, “os? cos39 + cos 4 usw. Ist. 


2 5 16 


» 
| 
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In a“ a, ausgedrückt lauten also nach der Bedeutung für die 1, die nicht verschwinden. 
bisi=4 


| a,” 
tan) + an” +12(1 
12 (1 Ay) | + Ay) +4 (6 Ay) A +3 Ay) 
Pr + a) 
12(2+a,’) (4 +a,)+4 (2 (4 +9 Ay") 
12(12-+a,) PR. 6-+a,) (12 + 
( Ay) (4 Ay + Ay) ( + Ay) ( | Ay)” 
| 4(6 + d+3ay?) 
2(16 4 Ay) 32 (1 + Ay) (A + (*) T au) (16 + 


Man erkennt, daß die Größen «a,,; eine für die a, konvergente Entwicklung sowohl in 
bezug auf » als auf k ergeben. 
Die Koeffizienten a;= SIa,,;y} erscheinen nun bis zur vierten Potenz von y in der 
foleenden Gestalt: 


für 4 bei Vernachlässigung der Potenzen von höher als +. 


Aus den a) folgen dann die Größen 5b, nach den Gl. (I1d) und zwar wird: 


- 


a” 


d, a. Ad, .a® 


für #>4 bei Vernachlässigung der Potenzen von z höher als 4. 
8. Die Druckfunktion. Es war in den Ansätzen (3) und (4) für die Verteilung des Druckes 
(genauer der mittleren Hauptspannung) gesetzt worden: 


wo nach (Ba), (7a) und (12) bis zu zweiten Potenzen von o: 


Die Einsetzung der oben erhaltenen Ansätze der a, liefert also zunächst: 
k 


n l 


1 + 
») 0° 


- 


+ sin (29) + + sin 0”) | 


16-+ a,” | 
> 


+sın + 


I 
0) | 
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end noch einfacher für 0= 0, d.h. am inneren Rande (p=-p;) 


(18a). 
+ sind a, 
‘ kann zunächst die beulen Grenzfälle sehr breiten und sehr schmalen Lagers betrachten. 
Der erstere Fall muß in der Umgebung von z=0, d.h. in der Nähe der Symmetrie- 
ne des Lagerzylinders auf den früher durehgearbeiteten, ebenen Fall (ZAMM 1935, S. S1, 
t mit Hilfe eines Grenzüberganges (a, >V) zurückführen, während für die Umgebung 
L. = n 
(1. | = Randquersehnittes ein in der Längsrichtung z andersartiger Lösungsansatz gemacht werden 
te. 


9, Zurückführung auf den Fall des schmalen Lagers. Der zweite Fall «,. 1(, >o) 
n 


ubt, die Ergebnisse der Formeln (15). (17) und (18a) erheblich zu vereinfachen. 


Dabei ist zu beachten, daß es in den Formeln (15) und (17) nur auf die Kleinheit von 
«, ankommt, die praktisch schon für kleine » erfüllt sein wird. 


ohl in Es war nämlich: 
| 
n za) (2n-+1)’ 
Verlangt man z. B., daß schon für »=2 bzw. =1, a,’ 100 sei, so wird die verhältnis- 
(16). | mäßige Lagerlänge 
0,7854 bzw. 0471. 
a 4 20 
Begnügt man sich mit a,’ =49 für n=2 bzw. =1, so ergibt sich: 
Da 3n 
== 1,122 bzw. —= 0.673. 
2a 14 11 


Es wird also im allgemeinen genügen, die Lösung für «, ?>V aufzustellen und sie mit Hilfe 
der Glieder für n=0 und »—=1 zu korrigieren. Häufig wird auch die Korrektion für n 0 
alleın genügen. 


(17). 


Für @,„ ?—>0 vereinfachen sieh nun die Formeln (15), (16) und (17) wie folet: 
\ 


.(lda), 
3 5 
.(16a), 
_ 5 
| ( (Ta). 


Mit den Werten für a, erhält man für die Druckfunktion am inneren Rande bei sehr 
nalem Lager und bis zur vierten Potenz der spezifischen Exzentrizität z 


- 


ww 


\ ( 


2 
\ 
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Die Summe der rechten Seite bedeutet nun eine parabolische Verteilung parallel 


.) 2 
/ylinderachse von der Form ll 2 mit dem Scheitel ın der Svmmetrieebene  ı 
I.) 

Ordinaten Null an den Rändern. 

Beachtet man noch nach S. 275, daß 9 gesetzt waren, so 


hält man in wirklichen Längen z,. d, und wirklichen Drucken p, für den Druck ps, ; 
inneren Rande: 


(1S])), 
7 d 2 


Die Korrekturglieder I pn,; werden, wenn man das Lager so schmal voraussetzt, daß 
man nur für korrigieren muß und wenn man die a;,;, für n=0 mit bezeichnet: 


Ta 


+) 


wo die Werte der a,;"” aus den Formeln (15) für n=-0 zu berechnen sind. 


10. Die Strömungsgeschwindigkeiten. Die Formeln ©), (4). (10) und (13) führen auf die 
folgenden Werte der tangentialen spezifischen Geschwindigkeit «, und der axialen spezifischen 
Geschwindigkeit (siehe Gl. 


„ k 
| k 


Diese länes der radialen Koordinate wechselnden Geschwindigkeiten sind der Beobachtung 
nieht zugänglich. Dagegen werden im wesentlichen, wie im nächsten Abschnitt gezeigt, «die 


SL Ya 
spezifischen Durchflußmengen (do und \wdo beobachtet. Diese sind: 
\ \ cos a, 2 | cos | + 
(20 
\ 
E (< + 2) Ay > 0 / 


In «diese Gleiehungen sind nun die Werte für o.. a; und b, aus den Gln. (19), (- 
(16) und (17) einzusetzen. Man erkennt, daß der allgemeine Fall explizit zu recht umstäı 
lichen Ergebnissen führt. 


Es möge deswegen hier nur dasjenige Ergebnis explizit hingeschrieben werden, welel > 
man erhält, wenn man sieh auf genügend schmales Lager und sehr kleine Exzentrizi 
beschränkt, d.h. also, wenn man die vereinfachten Formeln (15a), (16a) und (17a) verwen“ 
und in ihnen nur bis zur ersten Potenz von z geht. 


(1S]h). 


t, daß 
+ 


ruf die 
schen 


ehtung 


ot, (die 


(1: 
ıstäl 


releh 
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Die Gleichungen vereinfachen sich damit wie folet: 


va 

Ya 


M 
Die hier auftretenden trigonometrischen Summen stellen innerhalb der Lagerlänge I 
2 
Konstante 7 ‚ dann die Parabel — | schließlich die Grade 
d \ 
Geht man außerdem noch zu den wirklichen Geschwindigkeiten und (Gl. (1b)) 
| dem wirklichen radialen Abstand Ar= oa über und dividiert, um von Durchtlußmengen 
Geschwindigekeiten zu kommen, dureh die mittlere Liehtweite so ergibt sich schließlich 
für die mittleren Geschwindigkeiten und 


dar. 


477 
= ), \ 24608 | (l-+ycos)) 
| 
(20a). 
va 
ya 2, 


Bei der hier beschränkten Ordnung der Näherung gewinnt man schließlich die Gleichung 
ler Stromlinien aus: 
ds, > 


af 


und damit, wenn die 2-Koordinate einer Stromlinie bei = — bedeutet: 


zo e-t | 754% 608 . Pr . . . . . (21a). 


Die Stromlinien sind also Wellenkurven mit der Periode 27, deren 
\mplituden von der Symmetrieebene des Lagers nach außen von Null aus 
linear wachsen. 


Auch bei diesen Kurven kann die Näherung durch Korrektur auf größere Lagerbreite, 
nalog wie für die Druckfunktion, weiter getrieben werden, wobei jedoch die Integration der 
(21) entsprechenden Differentialgleichung schwieriger wird. 


11. Vorläufiger Bericht über eine Versuchsanordnung. Auf Anregung des Verfassers wurde 
in Institut von Herrn Kucharskıi ein Apparat folgender Anordnung konstruiert und gebaut: 
In einem sehr genau hergestellten Glasrohr von der liehten Weite von SO’ mm und der 
länze von 100 mm, das durch Querwände abgeschlossen ist, ist mit paralleler Achse ein 
Stahlzylinder von ganz wenig geringerem Durchmesser so gelagert, daß er quer zu seiner 
\chse etwas verschoben werden kann. 
Er ist etwas kürzer als der umschließende Glaszylinder, so daß an den Endwänden ein 
owisser Vorratsraum entsteht. 
Die Querverschiebung ist dureh zwei elastische Gummimembranen in den Querwänden 
nöglieht, welehe auch dafür sorgen, daß, nachdem der Zwischenraum zwischen den beiden 
/s\indern mit einer Schmierflüssiekeit ausgefüllt ist und an den Enden je ein Flüssigkeits- 
ratsraum entstanden ist, daselbst im wesentlichen der äußere Atmosphärendruck herrscht. 
Der innere Zylinder kann nun mit Hilfe einer Kurbel und eines Zahnradgetriebes in 
inelle Umdrehung versetzt werden, und wenn man der Schmierflüssigkeit Aluminiumpulver 
" Holzmehl zusetzt und durch die Wandung des Glaszylinders die Strömung beleuchtet, 
on man die mittleren Stromlinien, deren Bereehnung oben näherungsweise durchgeführt 
beobachten und photographieren. 
Der Apparat ist in Abb. 2 im Liehtbild dargestellt, und ferner ist ein Stromlinien- 
tbild in Abb. 3 wiedergegeben. 
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Es zeigte sich in der Tat, daß die Stromlinien von einfach pe... 
odischem Charakter nach den Stellen kleinen Druckes konvergieren u || 
von den Stellen großen Druckes divergieren, so daß sie also divergiere: | 
in die Seitenbehälter laufen und aus ihnen konvergierend wieder auf den 
Umfang des Zylinders erscheinen. 

Trotzdem die gewählte Anordnung kein schmales Lager darstellt, haben die Stromlin: n 
durchaus den im vorigen Abschnitt berechneten wellenförmigen Verlauf mit nach aubn 
wachsender Amplitude. Die unmittelbare Beobachtung zeigt diese wellenförmigen Str. 
linien besonders klar, während bisher die photographischen Aufnahmen infolge der Reflexion 
und der nieht genügend starken Kontraste nicht so deutlich sind. Übrigens hat das im Bil. 
sichtbare querverlaufende System einiger Linien nichts mit der Strömung zu tun, sondern si 
dureh Kratzer im Glas verursacht. 

Bei den Versuchen zeigte sieh bei starker Exzentrizität und Geschwindigkeit ein 
Luftblasenbildung, d. h. Zerreißen der Flüssigkeit teils dureh Unterdruck, teils durdı 
Kapillarität an der Stelle der stärksten Unterdrücke. 

Alle Versuche, auch die über Druckverteilung, Kavitationseintritt und Temperatur 
verteilung sind von Herrn Vogelpohl durchgeführt und von ihm auch weiter fortgeführt 
worden, worüber er in einer besonderen als Forschungsheft des V. D. I. erscheinenden 
Arbeit berichtet. Seiner Arbeit ist allerdings bisher nur die Reynolds-Michellsche 
Theorie der graden Keilströmung zugrunde gelegt. Es wird später auf die Folgerungen aus 
der hier vorgetragenen Theorie noch einzugehen sein. 


Abb. 


12. Schlußbemerkung. In einer späteren Mitteilung sollen die hier gegebenen Ansätre 
auch für teilgefüllte Ringbereiche durchgeführt und in Vergleich mit den Ergebnissen der 
ebenen Strömung von Mitt. 1’) und der räumlichen Keilströmung von Michell gesetzt werden. 


Verfasser ist sich bewußt, daß das vollständige Problem der Schmierung die Einflüs>« 
der ortsabhängigen Temperaturerhöhung durch Reibung, Konvektion und Leitung und «ir 
bei jedem Schmiermittel verschiedene Zähigkeits-Temperaturfunktion enthalten muß, jede 
schien auch sehon eine vollständige Durchführung der hier behandelten ersten Grundaufga'r 
qualitativ richtige Einblicke und Ansätze für die theoretische und experimentelle Erfassu: 
zu liefern, die zur Berücksiehtigung der oben genannten Einflüsse nun erweitert werd! 
müssen. 


’) In den Gl. (15a) und (15b) von Mitt. I sind drei Potenzexponenten der Inversionskonstanten k falsch angege! 
wie man bei Berücksichtigung der vorhergehenden Gleichung sofort einsieht. 
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'nfluß der Änderung einer Matrix auf die Lösung des zu- 
‚ehörigen Gleichungssystems, sowie aufdie charakteristischen 


Zahlen und die Eigenvektoren. 
Von Helmut Wittmeyer ın Seelze. 


Einleitung. Die nachfolgende Arbeit ist eine gekürzte Fassung der Dissertation des 
»rfassers mit dem gleichen Titel (Darmstadt 1934). Während in der Dissertation die Unter- 
chungen für komplexe Matrizen durchgeführt sind, werden in dieser Arbeit der Einfachheit 
‚ber sämtliche vorkommenden Matrizen durchweg als reell vorausgesetzt. Bei den 
‚oisten Sätzen sind die Beweise fortgelassen. Wo Beweise gebracht sind, sind sie in mehreren 

"ällen geometrisch gehalten — abweichend von der Dissertation. 

Die zuerst gefundene Abschätzung der Arbeit ist die Formel des Satzes I, $ 7. Sie 

»rmörlicht es, die Größe der Änderung der Lösung eines linearen Gle ichungssystems ab- 

‚uschätzen, wenn die Größe der Änderung der Koeffizienten bekannt ist. Auf die Bedeutung 

»iner solehen Abschätzung für die Bauingenieure hatte mich Herr Professor Dr. U. Wegener, 
Darmstadt, aufmerksam gemacht (siehe auch $ 11). Anschließend hieran ergaben sich dann 
init Hilfe ähnlicher Überlegungen die übrigen Abschätzungen. Auch sie haben Bedeutung für 
praktische Anwendungen: 

Satz IV, $S, z. B. ermöglicht eine Abschätzung der größtmöglichen Änderung der Eigen- 
(requenzen eines Systems von gekoppelten Massenpunkten bei Änderung der Kopplung zwise 
den sehwingenden Punkten. Satz V,$ 8, ermöglicht z. B. Fehlerabschätzungen bei der Be- 
reehnung von KEigenvektoren. 

In allen in dieser Arbeit gefundenen Abschätzungen kommen charakteristische Zahlen von 
vegebenen Matrizen (oder Matrizenprodukten) vor. Die Formeln wären natürlich praktisch 
unbrauehbar, wenn man in sie die Werte für die charakteristischen Zahlen wirklich einsetzen 
müßte. Nun kommen diese Zahlen in den Abschätzungen stets in der Weise vor, daß die 
\bschätzungen riehtig bleiben, wenn man in ihnen die absolut kleinsten charakteristischen 
Zahlen dureh absolut kleinere Werte und die absolut größten charakteristischen Zahlen durch 
absolut größere Werte ersetzt. Zur praktischen Anwendung dieser Formeln genügt also die 
Kenntnis geeigneter Abschätzungen für die charakteristischen Zahlen. 

Abschätzungen der absolut größten charakteristischen Zahl einer Matrix 
nach oben mit Hilfe ihrer Elemente sind schon in der mathematischen Literatur vorhanden 
's. Schur: Math. Ann. 66, S. 485). Weitere für unsere Zwecke besonders geeignete sind noch 
hinzugefügt worden (siehe $$ > u. 6). 

Abschätzungen derabsolutkleinstencharakteristischen Zahl einer Matrix 
nach unten dureh ihre Elemente sind in der Literatur noch nicht zu finden. Abschätzungen 
(dieser Art hat mir Herr Professor Dr. Wegner aus einer früheren, bisher nieht veröffentlichten 
Arbeit freundlicherweise zur Verfügung gestellt*). 

Bei der Untersuchung (reeller) niehtnormaler Matrizen A kommen in unseren 
Kormeln die größte und die kleinste eharakteristische Zahl des Matrizenproduktes 4’ 4 vor. 
\uch für diese Zahlen sind Abschätzungen durch die Elemente von .| gefunden worden. 

Im 3. Kapitel wird gezeigt, wie die Abschätzungen des 2. Kapitels mit Hilfe der Un- 
zleiehungen in $ 6 für numerische Rechnungen verwertbar sind. Dort sind zwei Beispiele 
‚ur Anwendung des Satzes I, $ 7, durchgereehnet. Die oben außerdem ee RR 
der Sätze IV und V sind in dieser Arbeit noch nicht durchgeführt worden * 


Bezeichnungen. 
kleine Buchstaben in Fraktur = reelle Vektoren im »-dimensionalen Raume. Ihre 
Komponenten werden mit den zugehörigen kleinen Buchstaben in Antiqua be- 


zeichnet, B. Zu). 
DB, (außer = große Buchstaben in Antiqua= reelle quadratische, »-reihige Matrizen. 


E ist die Einheitsmatrix. Die Elemente der Matrizen werden mit den zu- 
gehörigen kleinen lateinischen Buchstaben bezeichnet, z. DB. 


A (a; 


") Anmerkung bei der Korrektur: Wie mir Herr Professor Dr. Wegner mitteilt, ist eine seiner Abschätzungen, 


ich (1), 8 6, schon vorher von S. Gerse heorin: Bulletin de VAcademie des sciences de PUSSR (1051), S. 749, 
esen worden. 


") Anwendungen der Fehlerabschätzungen auf Iterationsverfahren s. S. 301 ff. dieses Heftes. 


j 
| | 
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1 zu A transponierte Matrix. Man erhält sie aus 4, indem man Spalten und Zeilen ı 
tauscht. 


4 Det. A Determinante, die als Elemente die Elemente der Matrix A enthält. 
Eine reelle normale Matrix ist eine solehe Matrix 4, für die eilt 4-44}. 
/u den reellen normalen Matrizen gehören 
I. die reellen symmetrischen, für die ja gilt A’—= A, 
2. die reellen schiefsymmetrischen, für die ja gilt A’ — A. 
/-! irgendeine eharakteristische Zahl einer Matrix A. Sie ist Lösung der „eharakie. 
ristischen Gleichung“ E!=0. Die charakteristischen Zahlen sind die 


proken Werte der „Eigenwerte* z, die Lösungen von |E- 2 A|=0 sind. 
7 ax Maximum der absoluten Beträge der charakteristischen Zahlen der Matrix I 
Entsprechende Bedeutung hat 
größte charakteristische Zahl der Matrix usw. 


Vektor mit den Komponenten: a, +... + 


. . . . 
Quadratische Form = Ihr wird die (symmetrische) Matrix A=(a; 
i,k | 
zugeordnet. Mit Hilfe dieser Matrix läßt sie sich abgekürzt schreiben: 
Ag; 
Die charakteristischen Zahlen 4, der «quadratischen Form F sind die 
charakteristischen Zahlen der zugehörigen Matrix A. Durch eine Haupt: 
achsentransformation geht die quadratische Form über in 


1,2, +42. +... + Ay - 


Kigenvektoren r einer Matrix .|, die zu einer charakteristischen Zahl 4 von A gehören, sind 
Lösungen des homogenen linearen Gleiehungssystems (4 44 E)r = 
mit der Nebenbedingung =1. 

Wiehtige Formeln sind in den Kapiteln durchlaufend mit römischen Ziffern numeriert. 

Die Numerierung mit arabischen Ziffern beginnt in jedem Paragraphen von vorn. 


I. Kapitel. Vorbereitende Abschätzungen. 
$ 1. Die Matrix .1’.] und ihre charakteristischen Zahlen '). 
I. Für beliebige reelle Matrizen A ist 1A eine reelle symmetrische Matrix. Ihr ist die 


positive reelle quadratische Form 


r-AAr=Ar- Are la. +... u) 

zugeordnet. Nach Durchführung einer Hauptachsentransformation erkennt man: Die Matrix 1 
hat für beliebige Matrizen „| nur nichtnegative charakteristische Zahlen. Hat A eine nichtver- 
schwindende Determinante, so sind die charakteristischen Zahlen von A’ A sogar sämtlich positiv. 

2. Bei geeigneter Numerierung der charakteristischen Zahlen einer reellen normalen 
Matrix .1 besteht zwischen ihnen und den charakteristischen Zahlen der Matrix AA die 
Beziehung 


3. Für reelle Matrizen 4 gilt?) 


14 1A 20 


I, Verl. z. B. ©. ©. Mae Duffee: The theory of Matrices, Ergebnisse der Mathematik und ihrer Grenzgebı!®, 
Bd. 2 (1035), oder R. Conurant und D. Hilbert: Methoden der mathematischen Physik, Bd. 1 (1931), Kapitel 1. 


®, E. T. Browne: Bulletin of the American Math. Soc. 34 (198), S. 363 bis 368. 


n 
Für reelle Matrizen ist gleich dem Maximum der Bilinearform unter den Ne! 
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$ 2. Die grundlegenden Ungleichungen. 


I. Aus 
pt 


>. Zu jeder der Ungleiehungen (I) läßt sich ein Vektor x* mit Ix*|=-0 angeben, für 
zusammen mit dem Vektor Ax” in der Ungleichung das Gleichheitszeichen steht. 


© Ungleichungen ergeben um so schärfere Abschätzungen, je dichter 4, und 4, zusammen- 
ren. 

>. Zum Beweise des in den Absätzen 1 und 2 Behaupteten, den wir — abweichend von 


Dissertation hier geometrisch im »-dimensionalen Raume durehführen wollen, beschränken 
» uns zunächst auf den Fall. daß die Matrix A eine niehtverschwindende Determinante hat. 
(1) folgt durch Quadrieren 


At 


ın ist allgemein Ux-3=xr- €" 3, wovon man sich leicht durch Ausreehnung überzeugt. Hiermit 
ergibt sich 


\uf der rechten Seite dieser Gleichung steht eine positiv definite, reelle quadratische Form 
in. den Komponenten des Vektors x mit der Koeffizientenmatrix 1A. Diese Gleichung bedeutet 
scometrisch: Zeiehnet man die zu allen y mit konstantem 'y gehörigen Vektoren x in ein 
orthogonales Koordinatensystem «derart ein, daß die Anfangspunkte der x-Vektoren in den 
Nullpunkt fallen, so liegen die Endpunkte dieser Vektoren auf einer Fläche 2. Grades. Diese 
"liche ist niehtausgzeartet und liegt ganz im Endlichen. (Im dreidimensionalen Raume hätten 
wir ein Ellipsoid vor uns.) Die Länge der x-Vektoren ist höchstens gleich der größten Haupt- 
achse dieser Fläche und mindestens gleich der kleinsten. Die Transformation der quadratischen 
"orm (2) auf Hauptachsen 


wo die 4; die charakteristischen Zahlen der Matrix A’A sind ergibt, daß die größte und 
die kleinste Hauptachse die Längen 


1 1 
haben. Daher gilt die Abschätzung 
*, 


die man als den Spezialfall der Ungleiehung (Il) für Matrizen A mit nichtverschwindender 
Determinante von A erkennt. 

Für einen Vektor x”, der die Riehtung der kleinsten Hauptachse hat, und den nach (1) 
Ilm zugeordneten Vektor 9” gilt in dem linken Teile der Ungleichung offenbar das Gleich- 
oitszeichen. Entsprechendes gilt für einen Vektor x*, der die Riehtung der größten Hauptachse 
hat. Die eben angeführten Vektoren x* sind, wie man leicht erkennt, solche Vektoren, die 
len Eigenvektoren der Matrix A’A proportional sind, die zu der größten und der kleinsten 
ohurakteristischen Zahl dieser Matrix gehören. Diese Vektoren x* lassen sich daher stets so 

bringt man die oben bewiesenen Ungleiechungen auf die Form (l), so kann man zeigen, 
\ub aus Stetigkeitsgründen das in den Absätzen 1 und 2 Behauptete auch noch für Matrizen 

 versehwindender Determinante gilt. 


S 3. Schranken für die Absolutwerte der charakteristischen Zahlen der Matrix 1: BR. 


l. Abschätzungen nach oben. 


a) Es ıst?’) 
„4A „BB 
/ / . 


"max "max "max 


') Anmerkung bei der Korrektur: Die Abschätzung (ID fand ich nachträglieh in etwas anderer Form aueh bei 
ert: Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen, Leipzig 1912, S. 120, und bei 
Intner: Spektraltheorie der unendlichen Matrizen, Leipzig 1929, S. 150, wo auch die Abschätzung (VIID an- 
ırt 
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Beweis: Bei der Überführung des Vektors x in den Vektor y dureh die Transformation v- 
wird der Vektor x zufolge der Ungleiehung (I) höchstens auf das Full -fache gestre: 
Ähnliches gilt für die Transformation 3 Ay. Werden die beiden genannten Transformatioı 
nacheinander ausgeführt, so geht dadurch der Vektor x in den Vektor 3 ül er und wird höchst«: - 
auf das -fache gestreckt: 


"mi 


Nun können wir die Transformationen B und A in die eine Transformation 4 B zusammenfass: 


Nachg 2, Absatz 2, gibt es einen Vektor x” mit x”) ->0 und einen zugehörigen Vektor "= AB ı 
für die gilt 


"MAX 


Setzen wir die Gl. (3) in die Gl. (1) ein, die ja auch für 3= 3" und x = x* richtig ist, so or 
halten wir 


Nach Division dureh x” >00 erhalten wir hieraus die behauptete Ungleichung. 
b) Es ist er 


Beweis. Nach $ 1, Absatz 5, ıst 


MAX "max 
Wird dies in (Ib) berücksichtigt, so ergibt sieh (III). 
e) Für (reelle) normale Matrizen A und B gilt 


max Max max 
Beweis. Man berücksichtigt in (III das in $ 1, Absatz 2, Gesagte, 
>» Absehätzuneen nach unten. 


Die Beweise der folgenden Ungleiehungen sind fortgelassen, da sie denen des Absatzes | 
eanz ähnlich sind. 


»BB 
mm min 
mim 


"ir (reelle) normale Matrizen I und B gilt 


‚IB .I2B 


nın 


$s 4. Schranken für die Absolutwerte der charakteristischen Zahlen der Matrix (.1+B). 


Der Beweis der folgenden Formeln kann ebenso wie in $ 3 durch einfache geometrische 
Überlegungen erfolgen. 


I. Abschätzungen nach oben. 


max 
Max 


ür (reelle) normale Matrizen und eilt 


jA+B < 124 ;B \ 
2. Abschätzungen nach unten. 
max 
min 


Für (reelle) normale Matrizen A und B gilt 


— [AB . (Al 


IM. tober 1936 Wittmever, Einfluß der Änderung einer Matrix 9] 


Bemerkung. Daß die Ungleichungen (IV), (VII), (X) und (XTIl) für reelle niehtnormale 
trizen nicht allgemein riehtig sind, erkennt man leicht an Gegenbeispielen. Die Ungleichung (IV) 


io1 B. nieht richtig für 
Stens 01 0 
dann ist 
ISS 


- 5. Abschätzungen der charakteristischen Zahlen der Matrizen .| und .1’.I durch die Elemente 
IBı*, | der Matrix 
I. Es ist 
4 
1 


SO 


. 
mim 


Beweis’): 4 sei irgendeine charakteristische Zahl von A, also irgendeine Lösung von 
Dann hat das Gleichungssystem mit der Matrix (A EA): 


n 
— a); — 00, Gz=l,...,) 
k 
niehttriviale Lösungen. Es sei 
e=1,...;#. 
Dann wird 
N 
k—1 
ki . 
und daraus folgt 
n, 
Ilieraus ergibt sieh sofort die oben behauptete Ungleichung. Sie hat natürlich nur Bedeutung, 
m wenn die rechte Seite größer als Null ist. 
| 2, Es ıst 
n 
(VIn. Beweis: Aus der Ungleichung (2) folgt 
n 
B). A| ii 
trische 
und schließlich 
| | 
/ +2 


"u z. irgendeine eharakteristische Zahl von A ist und © ein dazu passender Index, so können 
(IX). rin dieser Ungleichung beide Seiten dureh ihr Maximum ersetzen, womit die Ungleichung (3) 
bewiesen Ist. 


». Weitere Abschätzungen finden wir mittels der Ungleichungen (VIIT) und (Ah) aus $ 4: 


(A) - 
(N a) Wir zerlegen die gegebene Matrix A in folgender Weise: 
(XI! », U, Wegner: siehe Einleitung und Anmerkung *). 


1y* 


(\ 
| 


Zisehr. f. ang 
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wobei 


— Ay; 


IV 


ist. S und 77 sind als reelle symmetrische bzw. schiefsymmetrische Matrizen normal. \ 
setzen nun in den Ungleiehungen (4) und (5) 


B=S und C=UD 


Da S und U normal sind, können wir das ins$ 1, Absatz 2, Gesagte anwenden und erhalten 


(; 
\"min 


. WL 1728] 
4 min max \ max 4 max | max * 


Zur weiteren Abschätzung benutzen wir die aus (1) und (3) folgenden Ungleichungen 


n n 
Min — 2 AU| max Max PAUSE 
Dann ergibt sich 
(A Min (2 = |a;k + — Max la: . (6), 
n n 


P) Weitere Abschätzungen erhalten wir dureh folgende Zerlegung der Matrix 4: 


A=B-+ÜC, 
Wo 


— für i=k für ick 


| 
0 für | 


h; . 
| für 


Da die Matrix PB eine Diagonalmatrıx ist, gilt, wie man leicht einsieht, 


Max "min 


Weiterhin gilt die Abschätzung 


n 


die leieht aus einer bekannten Abschätzung von Schur’) folgt. 


Setzen wir nun in (4) und (#9) B+ CA und berücksichtigen (S) und (9), so erhalten wir 


l 


$ 6. Zusammenstellung der Abschätzungen des 5. Paragraphen '‘). 


Bei der folgenden Aufzählung der im letzten Paragraphen bewiesenen Abschätzung«' 
wird dort, wo es praktisch wichtig erscheint, gleich die Ungleichung 


4°“ \min AA Max 
aus $ 1, Absatz 3, berücksichtigt. Die in den folgenden Formeln enthaltenen Abschätzung ' 
für die charakteristischen Zahlen der Matrix A4’.4 sind durchweg neu. Die beigefügten Zitı\" 


eben an, von wem die in den Formeln enthaltenen Abschätzungen für die charakteristisch 
Zahlen der Matrix .| zuerst angegeben sind. 


4, ]). Schur: Math. Ann. 66 (1909), S. 492, 2. 
4a) Weitere Abschätzungen folgen in einer späteren Arbeit. 
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1. Abschätzungen nach unten 


4 
l k l 
3, Abschätzungen nach oben (i=|],....., n). 
i-Ek 
2 Max (2 + Max (2 Ur — Aki (6). 
Ü l 7 
2. Kapitel. Die Hauptsätze. 
$S 7. Die Lösungen benachbarter linearer Gleichungssysteme. 
Satz I. Es seien die beiden linearen Gleichungssysteme gegeben 


ER 


Die (reellen) Matrizen A und B mögen der Ungleichung 


"max/ 


senügen. Dann gilt folgende Abschätzung für den „Unterschied* der Lösungen x und x* der 
„benachbarten“ Gleiehungssysteme (1) und (2): 


1 
(ax) 


Beweis. Es ist 
d)=(A+B)r"— Ar—Br. 
Die Beachtung der Gl. (1) und (2) ergibt 
Bı=f-—B:. 
Schätzen wir mit Hilfe der Ungleichungen (D, $2, die linke Seite dieser Gleichung nach 


unten ab und die rechte Seite nach oben, so erhalten wir 


"min 


ir beachten nun, daß ebenfalls mit Hilfe der Ungleiehungen (N, $2, aus (1) folgt 


"min 
ıterhin berücksichtigen wir die Ungleichung (XI), $ 4. So erhalten wir aus (4) 
ı 
ı «ividieren wir durch den Faktor von 'r  x* auf der linken Seite dieser Ungleichung, 
wegen (3) positiv Ist, und erhalten dadurch die zu beweisende Ungleichung (N). 


°) Bei der Benutzung der Abschätzungen (6) und (6) ist es zweekmäßig, sieh zunächst die beiden symmetrisehben 
) und ( ) aufzuschreiben, 


| 
l. Wir 
| 

| 
| 
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Bemerkune. Es fällt auf, daß man das Niechtverschwinden der Determinanten ı 
A und 4 -+ B nieht besonders vorauszusetzen braucht. Der Grund ist, daß dies schon ı 
selbst aus der Voraussetzung (3) folgt. Man erkennt dies leieht, wenn man die | 
eleichunge (XII), $ 4, beachtet. 


Zusatz I, I. Unter den Voraussetzungen des Satzes I gilt folgende Abschätzung | 
den „relativen Unterschied“ der Lösungen x und x” der „benachbarten“ Gleichungssysteme 
und (2): 


(„B’ 


min 


Beweis. Nach $ 2, Absatz 1. folgt ans der Gl. (1) die Ungleichung 
A\—} 

Durch Division der Ungleichung (1) dureh (5) erhalten wir (ID. 


Bemerkung. Die Abschätzungen (D und (ID) sind um so schärfer, je diehter 

und 


A 


zusammenliegen. Dies folgt daraus, daß die zur Ableitung dieser 


schätzungen benutzten Ungleiehungen (1), $ 2, sich ebenso verhalten. 


Zusatz I, 2. Sind die Matrizen .| und B normal, z. B. symmetrisch, so vereinfachen 
sich die Ungleiehungen (D und (ID zu 


;B\ 
£ 
Y min (N 
| 44 ‚B 
min MAX 
| ‚B 
4° 
< | 4“ min (IT) 


Ss. Vergleich der zu benachbarten Matrizen gehörigen quadratischen Formen, 
charakteristischen Zahlen und Eigenvektoren. 


Die zum Teil recht langwierigen Beweise der folgenden Sätze sind zwecks Raum- 
ersparnis fortgelassen. 


Satz 11. 


n n 
i,k l | 
und 


seien zwei reelle quadratische Ausdrücke. Die Vektoren 
seien also reell, und die Matrizen 1 = {a;,) und B= (b;,) reell symmetrisch. Ferner seı 


44 AB 


min max 


Dann eilt für die Differenz m — m*\ der stationären Werte m und m” dieser „benachbarten 
quadratischen Ausdrücke Fund F* die Abschätzung 


2 
Ir 2 


Satz III. 4 sei eine (reelle) normale Matrix, B eine beliebige reelle. Dann eilt 


"MAX 


.angı 
nd Me 


\ 


lichter 
Ah- 


fachen 


Raum- 


se] 
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Andere Fassung: In der komplexen Zahlenebene ist der Radius des kleinsten (größten) 

jses um den Nullpunkt, der sämtliche (keine der) eharakteristischen Zahlen der Matrizen- 
.. .. 1 . . 

ıne (4 + B) enthält, höchstens um Freu, ° erößer (kleiner) als der Radius des kleinsten 


sten) konzentrischen Kreises, der sämtliche (keine der) eharakteristischen Zahlen der 
‚malen Matrix 4 enthält. 


Satz IV. Es seien A und C=(A-+ B) normale Matrizen. Bei geeigneter Numerierung 

eharakteristischen Zahlen A,-1, der Matrix A und 4,©,...,4,° der Matrix € 
dann 


"max 


Andere Fassung: Sind A und © (.4-+B) normale Matrizen, so liegen in der komplexen 
'hlenebere die eharakteristischen Zahlen von € in Kreisen mit dem Radius Senf ae die die 
-lırakteristischen Zahlen von A als Mittelpunkt haben. Und zwar liegt in jedem Kreise 
sıindestens eine charakteristische Zahl von €. 


Satz V. A sei eine normale Matrix, 4 sei eine zu A gehörige, ein- oder mehrfache 
charakteristische Zahl. B sei eine beliebige reelle Matrix, deren Elemente jedoch absolut ge- 
nommen so klein sein sollen, daß es eine zur Matrix (A-+ B) gehörige charakteristische 
/ahl 4° gibt mit der Eigenschaft 


/ 
MAX 
Ilierin sollen die /; die von / verschiedenen ceharakteristischen Zahlen der Matrix A sein 


Dann gibt es zu jedem zur charakteristischen Zahl 4* gehörigen Eigenvektor x" der 
\Matrix (.1 + BD) einen geeigneten Eigenvektor x der Matrix A, der zu der charakteristischen 
Zahl 4 gehört, sodaß gilt 


Min — A| 
Diese Abschätzung kann auch durch die handlichere, jedoch weniger scharfe Abschätzung 
(VIII) ersetzt werden: 


| 


Mın A| 


Bemerkungen: 


I. Die Bedingung (1) schränkt die Anwendungsmögliehkeit der Formel (VII) nicht be- 
sonders ein, da, wenn die Bedingung (1) nicht erfüllt ist, die Formel entweder etwas Triviales 
oıler etwas Unsinniges aussagen würde. 

2. Die Unbestimmtheit, die in der Aussage: „Dann gibt es“ liegt, ist dadurch bedingt, 
(laß 2 eine mehrfache charakteristische Zahl sein kann. Wenn dies nämlich der Fall ıst, 
so gehört zu / eine Schar von Eigenvektoren. Wie man geometrisch leicht erkennt, erfüllt 
dann tatsächlich nicht jeder dieser Vektoren x unsere Ungleichung (VI). 


Wir vermerken noch folgenden Sonderfall des Satzes V: 


Zusatz V.1. Aund A® (14-+ B) seien zwei (reelle) normale Matrizen. A habe lauter 
einfache eharakteristische Zahlen, und es sei 


.) ®<d, wo d=Min 


Dann lassen sieh die ceharakteristischen Zahlen bzw. und «die zugehörigen 
„envektoren x; bzw. x;* der Matrix A bzw. 4* so numerieren, daß gilt: 


w a 


“B’ B\'la 


er einfacher: 


2y2 


e | 
(»)) 
in 
\ | | 
( 
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Bemerkungen: 
I. Ist .I eine Diagonalmatrix, so Ist 


d Min 


>». Ein einfach anzuwendendes hinreichendes Kriterium dafür, daß eine Matrix A= («a 
keine mehrfachen charakteristischen Zahlen hat, teilte mir Herr Professor Dr. Wegener n 
‚Die Matrix 4 hat nur einfache charakteristische Zahlen, wenn für alle © und k 


„ 
la; — >2 Max 2 [aus 
ist.“ 
3. Die Voraussetzung (2) sorgt übrigens dafür wie aus Satz IV ersichtlich —, daß 


auch die eharakteristischen Zahlen der Matrix A’= A-+B einfach sind. 


I. Die Ungleiehungen des Zusatzes V,1 ermöglichen es z. B., abzuschätzen, welchen 
Winkel einander entsprechende Hauptachsen zweier benachbarter quadratischer Flächen 
höchstens miteinander bilden, und um wieviel sich diese Hauptachsen höchstens in ihrer 
Länge unterscheiden. 


3. Kapitel. Anwendungen und Beispiele. 


Von den mannigfachen Anwendungen, tie die Abschätzungen des 2. Kapitels zulassen, 
seien nur einige herausgegriffen. 


$ 9. Grundsätzliches über die praktische Anwendung der Abschätzungen des 2. Kapitels. 


! | ı 9 


I. Die Abschätzungen für die Größen , m” |, des 2. Kapitels wird 
man in der Praxis in folgender Weise benutzen können: Die nieht mit einem Stern versehenen 
Größen wird man als wahre Lösung des ın Frage kommenden Problems auffassen und die 
zesternten Größen als Näherungslösung dieses Problems. Ist dabei die Näherungslösung 
wirklich als exakte Lösung eines „benachbarten“ Problems gewonnen, so sind unsere Ab- 
sehätzungen unmittelbar anwendbar ($ 10, Beispiel 2). Ist dagegen die Näherungslösung 
irgendwie in anderer Weise gefunden worden, so wird sich oft leicht zu dieser Lösung ein 
zugehöriges Problem (z. B. Gleiehungssystem) aufstellen lassen, dessen exakte Lösung sie ist 
($ 10, Beispiel 1). 

2. Bei praktischen Fehlerabschätzungen wird man natürlich die in den Ungleichungen 
des 2. Kapitels vorkommenden eharakteristischen Zahlen durch ihre in $ 6 angegebenen Al- 
sehätzuneen ersetzen. Daß sie dadureh richtig bleiben, ist schon in der Einleitung bemerkt 
worden. 


3. Kennt man von einem Problem mit der wahren Lösung = (ir, .#.......,) eine 


Näherungslösung (r,*, so kann man z.B. folgende Größen als Fehler der 
Näherungslösung x” bezeichnen: 


a) die in dieser Arbeit abgeschätzte Differenz 


b) «den maximalen Fehler in den Komponenten: 


Gel...) 


für ihn gilt, wie man leicht sieht, 
öo<4; 
e) den „durehsehnittlichen Fehler in den einzelnen Komponenten“: 


n 


\' 


“, In der Dissertation ist in dieser Ungleiehung versehentlich der Exponent , fortgelassen. 


| 
i—=] 


chr. f. ang 
th. und M: 


A=la 


enern 


eh daß 


welchen 
r Flächen 
s ın ihrer 


zulassen, 


\apitels. 


itels wir 
rersehenen 
ı und die 
ıneslösung 
ıinsere Abh- 
ıneslösung 
‚Ösung ein 


Ing sie ist 


leichungen 
benen Ab- 
bemerkt 


eine 
‘ehler der 
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ihn gilt nach der Schwarzschen Ungleichung 


| 
In 
d) den „relativen Fehler“: 


e) den „prozentualen Fehler“: 


— 


9,00 


ie Definitionen unter a), d) und e) dürften — nieht nur auf Grund der vorliegenden Arbeit 
ie für die Lösungen von linearen Gleichungssystemen angemessenen sein. 


$ 10. Zwei Beispiele zur Anwendung des Satzes I ($ 7). 


Erstes Beispiel. 
\ufgabe: 


ci eine gegebene Näherungslösung des (symmetrischen) Gleichungssystems mit überwiegender 
Hauptdiagonale: 


26,7 + 130,— 1,888 
S,lır, 620,+ 12x, = — 58,203 
130,4 620,4 — 430,— 88,390 (N). 
032,+ 12%, 430, 427,10, 98,664 
Wie groß ist die Genauigkeit dieser Näherungslösung”? 
a) Lösung: Wir schreiben die Gl. (1) in abgekürzter Form: 


Setzen wir die Näherungslösung x* in die linke Seite der Gleichung (1) oder (1”) ein, so er- 
eibt sich nieht Ay®=r, sondern 


Wo 
(0,598: 0.190: 0.764) 


ist. Die Differenz |'r — x*! können wir, da 4 symmetrisch ist, mit Hilfe der Ungleichung (1), 
$ 7, abschätzen. Hierin haben wir noch |4/P 0 zu setzen. So erhalten wir 


min 


(3). 
/" ,;„ können wir in dieser Ungleichung zufolge (1). $ 6, noch ersetzen durch 
Min (la;; — 

Oheses Minimum berechnet sieh aus den Koeffizienten der zweiten Zeile der Matrix A zu 
238,6 - 5,1+62+12)= 16,1. 
Veiter ist 

(0,538? 4 0,0032 + 0,190° 0,76P)"2 0,95. 


Ungleiehung (3) ergibt dann 


0,05 
"aus erhalten wir nach $ 9, Abs. 3, 
0,06 
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Um den relativen Fehler zu berechnen, benutzen wir die Ungleichung (I1’), $ 7, die ja 


symmetrische Gleichungssysteme Gültigkeit hat, und setzen in ihr wieder \4P x =0: 


44 
Nun ist nach 6, 
/A| Max 2 
| 


Das Maximum auf der rechten Seite bereehnet sieh in unserem Falle aus den Element: 
der zweiten Zeile von A zu 


51 +8B6+62 +12 


Ferner ist 
(13,888 58,203? 88,390? 4 98,66. P) 151. 


Nunmehr können wir die Ungleichung (7) zahlenmäßig auswerten und erhalten 


0,95 


une daraus 
. . . . . . . . . . . . . . . (9), 


100 
b) Untersuchung der Schärfe dieser Abschätzungen. 
Die exakte Lösung des Gleiehungssystems ist 
— (2,020; 2,500: 3,010; 3,970). 
Mit Hilfe dieser Werte berechnen wir die genauen Werte der Größen 4,6, 0,0, 9,0, ZU 


A= (0,02? +0,01? +0,03) =0,08 . . 2 2 2 2.2.4), 


0.06 
4 
“ 5.9 
= 0,65% (9. 


Der Vergleich der wahren Größe der Fehler mit ihren Abschätzungen ergibt, daß unsere A\b- 
schätzungen die Größenordnung des Fehlers gut getroffen haben. 


Zweites Beispiel. 


Aufzabe: Um für das (unsymmetrische) Gleichungssystem mit überwiegender Hauptdiagonale 


+0,05 — 001 l | 


+00 10) 
eine Näherungslösung x* (,*, .r,*) der wahren Lösung x = zu bekommen, 


man «das vereinfachte Gleichungssvstem 


0,05, | 
-0.03 - ( 


auflösen. das aus (a) dureh Nullsetzen der beiden kleinsten Koeffizienten der Unbekann! 
Welehen Fehler wird man höchstens begehen, wenn man die Lösung von 


entstanden Ist. 
als Näherungslösune von (a) ansieht” 
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ie ja | a) Lösung: 

IE 


setzen abkürzend für (a) 


| für (b) 


wird 


ementen 00 
nd 
/ur Abschätzung von - x#! benutzen wir die Formel (D, $ 7%, in der wir =0 setzen 
können: 
ie ( ) 
(S) (Amin } min 3 
ıst 
\ach Formel (5), $ 6: 
H 
i=1,..+9 i,k=1 
ik 
ist 
(N), Nach Formel (2), $ 6: 
i,k=]1 
i-Ek 
(6), 
Wo 
(0,05? 0,012 4 +0,08 + 0,09 (h) 


ist, also 
sere Ab- 


min} 
\\ır setzen (fJ) und (i) in (e) ein und erhalten 


agonale 109.3:0.0142 


Oamaus ergibt sieh weiter (vgl. $ 9, Absatz 3): 


0,19 
o< 
/ur Abschätzung des relativen Fehlers benutzen wir die Ungleichung (ID), $ 7, in der 
" gleich setzen: 
IX A: 44 2 12 
16 
Auswertung dieser Ungleichung haben wir nur noch die Größe "als '? abzuschätzen, wo- 
ann! wir die Formel (e) benutzen. Bei Beachtung von (y) und (h) ergibt sieh 
von 
u=1-+ 009 = 1,094 


0) 0 
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Wir setzen nun die Werte aus (f), (D und (i) in (k) ein: 


1.09 


. -< 002. 
00015 


Hieraus folet 


b) Untersuchung der Schärfe dieser Abschätzungen. 


Die Lösuneen von (a) und (b) sind 


r = (0,804; 1,925; 9,919), 
bzw. 
= (0,899; 2,027; 9,914). 


Hiermit berechnen wir die wahren Werte von J, d,6,0, und vergleichen sie mit den 
oben erhaltenen Abschätzungen. 


Wahre Werte: 


Abschätzungen: 


100 


Wir sehen wieder, daß die Fehler der Größenordnung nach gut abgeschätzt sind. 


$ 11. Vergleich der Fehlerabschätzungen des 2. Kapitels mit früheren Fehlerberechnungen. 


Fehlerbereehnungen für Lösungen linearer Gleichungssysteme, und zwar für solche, die 
bei der Berechnung mehrfach statisch unbestimmter Systeme vorkommen, sind nach einer 
eanz anderen Methode schon von Hertwig’) und Pıirlet°) durchgeführt worden. Die Resultate 
unterscheiden sieh von den unsrigen im wesentlichen in folgenden Punkten: 


I. In den dortigen Arbeiten wird der Fehler (der als sehr klein vorausgesetzt ist) an- 
eenähert berechnet, wobei Glieder von höherer als erster Ordnung vernachlässigt werden. 
Da jedoch die Koeffizienten der Glieder von höherer als erster Ordnung sehr groß werden 
können, bietet diese „angenäherte Berechnung“ des Fehlers keine wirklichen Anhaltspunkte 
über seine wahre Größe. 


Der Satz I, $ 7, bietet dagegen die Möglichkeit zu sagen, wie groß der Fehler höchstens 
ist. Über seine Größenordnung braucht dabei nichts vorausgesetzt zu werden. 


2. Die „angenäherte Berechnung” des Fehlers nach der Methode von Hertwig gestaltet 
sieh im allgemeinen sehr schwierig. Hertwig selbst sagt darüber (l. e., S. 195, 2. 61.): 
“.... Jedoch wird die Bereehnung noch zeitraubender wie die Lösung der Gleichungen.“ 
Man muß nämlich dabei die Svstemdeterminante und sämtliche Unterdeterminanten (n - Liter 
Ordnung kennen. Die Auswertung dieser Determinanten gestaltet sich dann einfacher, wenn 
die zu bereehnenden Fachwerke viele Symmetrieachsen haben. 


Unsere Fehlerabschätzung ist im Falle eines allgemeinen Gleichungssystems zwar auch 
etwas schwer zu handhaben, da man dann zur Abschätzung der in ihr vorkommenden eharak- 
teristischen Zahlen auch eine Abschätzung®) benutzen muß, die die Berechnung der System 
determinante voraussetzt. Immerhin ist die Berechnung auch hier noch viel einfacher als 
die Bereehnung in den oben erwähnten Arbeiten. In dem praktisch wichtigen Falle, in den 
das Gleichungssystem eine Koeftizientenmatrix mit überwiegender Hauptdiagonale hat, ist «1 
Fehler sogar sehr leicht abzuschätzen, wie wir in $ 10 gesehen haben. 


Auf das Fehlen einer praktisch brauchbaren, einwandfreien Fehlerabsehätzung für 0 
Lösungen linearer Gleichungssysteme, besonders für die Zwecke der Bauingenieure, macht 
mich Herr Prof. Dr. Wegner, Darmstadt, aufmerksam und regte so die vorliegende Arbeit on. 

') A. Hertwige: Festsehrift Adolph Wüllner, 1905, S. 194 ff. 


s), Joseph Pirlet: Fehleruntersuehungen bei der Berechnung mehrfach statisch unbestimmter Syst 
Dissertation Aachen 109. 


», Diese von U. Wegner herrührende Absehätzung findet sieh in der Dissertation auf Seite 35 als Forme! 
Sie ist in dieser Arbeit zwecks Raumersparnis fortgelassen. 
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Über die Lösung von linearen Gleichungssystemen 
durch Iteration. 
Von Helmut Wittmeyer ın Seelze. 
Einleitung. 


1. Ziel. In der Arbeit „Praktische Verfahren der Gleiehungsauflösung* von R. v. Mises 
„lH. Pollaezek-Geiringer in der Zeitschrift für angewandte Mathematik und Mechanik, 
9 (1929), S. 58 bis 77, untersuchen die genannten Verfasser folgendes Iterationsverfahren: 


+ Ca (x) 


sur Auflösung des linearen, inhomogenen Gleichungssystenis 


Sie stellen einige Konvergenzbedingungen auf und geben eine Abschätzung des Fehlers beim 
\bbreehen des Iterationsverfahrens beim “ten Schritt (Sätze I bis 5). Im folgenden sollen 
Iterationsverfahren für lineare Gleiehungssysteme von einem anderen allgemeinen Gesichts- 
punkte aus betrachtet werden, wodurch einige neue Ergebnisse erzielt werden. 
>, Bezeichnungen. 
1,B.0,... = große Buchstaben in Antiqua == reelle') quadratische, »-reihige Matrizen. 
Die Elemente der Matrizen werden mit den zugehörigen kleinen lateinischen 
Buchstaben bezeichnet, z. B. 


A=(a;}). 


E Kinheitsmatrix, 
4’ = die zu A transponierte Matrix, 
j.l irgendeine charakteristische Zahl der Matrix 4. Diese Zahl ist Lösung der 


„eharakteristischen Gleiehung*: Det. (4-4 


Se größte der (sämtlich niehtnegativen) charakteristischen Zahlen der Matrix A’A, 
Maximum der absoluten Beträge der charakteristischen Zahlen der Matrix 4. 


>. Wichtige Formeln’). 
a) Aus y= Ar folgt: 


"max \ min 16 
b) Für (reelle) normale (speziell symmetrische) Matrizen gilt 


1234| 


"max 
e) Es gelten folgende Abschätzungen für charakterıistische Zahlen ’) °*): 


| AA (A in (Min (2 k . . . . . (3), 
- IN 


') Die Behandlung komplexer Matrizen bringt keine weiteren Schwierigkeiten. 

°) Beweise finden sieh in: H. Wittmever: Einfluß der Änderung einer Matrix auf die Lösung des zugehörigen 
Vlehnngssystems, sowie auf die charakteristischen Zahlen und die Eigenvektoren, Dissertation Darmstadt, 1934 (im 
wenden zitiert mit „Dissertation*). Eine gekürzte Fassung der Dissertation findet sieh unter dem gleichen Titel 

IT. dieses Heftes. 

*) siehe Dissertation, Seite 6. 

%, siehe Dissertation, Seite 4. 

°) siehe Dissertation, Seite 35 bis 36. Dort finden sieh weitere Absehätzungen, auch solehe von A. Hirseh, 


chur und U. Wegner. Die in (3) enthaltene Absehätzung für „A min Ist zuerst von U. Wegner bewiesen 
den. Die in den Ungleiehungen (3), (4), (5) und (8) enthaltenen Ungleiehungen zwischen den eharakteristischen Zahlen 
Matrizen A und „PA rühren her von E. T. Browne: Bulletin of the American Math. Soe. 34 (1928), Seite 363 


. 


°a) Weitere Absehätzungen folgen in einer späteren Arbeit. 


Ztsehr. f.an: 
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n 
» (2Ja;;|— 2 + a, ;|)— Max 1}. 
i | i 1 
i 
n 
i | i 
[3 4 
AA |nax Max 21a} 
> Max 7\ 
IR l, ! | 
n 
i i,k 


$S 1. Aufstellung von hinreichenden Konvergenzbedingungen für ein Iterationsverfahren von 
allgemeiner Form. 


I. Das oben angeführte Gleichungssystem lautet in Matrizenschreibweise: 
mit 
und das Iterationsverfahren: 


In der oben zitierten Arbeit von R. v. Mises und H. Pollaczek-Geiringer ist in Gl. (10) 
für D eine Diagonalmatrix mit den Diagonalelementen e,,e,,...,c, gewählt. 
2. Wir wollen zeigen, daß das Iterationsverfahren (10) auch dann zum Ziele führt, wenn 
D eine beliebige Matrix — d. h. nicht unbedingt eine Diagonalmatrix — ist, sofern nur die 
Elemente von D gewisse „Konvergenzbedingungen* erfüllen. 
a) Wir setzen 
mit Determinante D==0, 
Dann ergibt das Iterationsverfahren (10) 


(E+- DA)” Dr—- (E+DA)Dr— Dr, 


(E+-DAliDr— Dr — Dr. . 
Falls die für © > entstehende Reihe konvergiert, wird 


r=limıd =— . . 


i>« 
Für dies x gilt dann: 
Y Dt, 
Dr, 
54:2 


D. h.: Konvergiert die Reihe (12) für irgendeine Matrix D mit nichtverschwindender Dei 
minante, so ist der Limesvektor x des Iterationsverfahrens (10) die Lösung des Gleichun > 
systems 


6) Diese Abschätzung ist in der Dissertation noeh nicht bewiesen. Sie läßt sich jedoch leicht aus der vo! 
gehenden folgern, wenn man beachtet, daß die Matrizen „# und „1’ dieselben eharakteristischen Zahlen haben. 

6a) Die unter Voraussetzung dieser speziellen Anfangswerte abgeleiteten Sätze I und II gelten ebenso für 
liebige Anfangswerte x'®, Die spezielle Wahl von x" = Dr vereinfacht nur unsere Beweise. 
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b) Wann konvergiert die Reihe (12)? 


Die Reihe 


vergiert bekanntlich, wenn < ist, die Reihe 
| 


ı 
E+B+B+B’+... 
dann, wenn 


| | 
> 
oder 


<1 ıst. 


Dies Kriterium wenden wir auf die Reihe (12) an: 
Die Reihe (12) konvergiert, wenn \ZAF+PA ax <t Ist, also wenn alle charakteristischen 
n von /alılen von (E+ D A) ganz im Einheitskreise der komplexen Zahlenebene liegen. 
e) Wie drückt sieh dies Kriterium durch die charakteristischen Zahlen von D A aus? 
Die charakteristischen Zahlen 4 von (E-+-D 4A) genügen der Gleichung 


Det. oder Det. (DA—- A 


(9), 


Die charakteristischen Zahlen « von DA genügen der Gleiehung 


/;wischen diesen Zahlen s und 4 besteht also bei geeigneter Zuordnung die Beziehung 

. (10), 

Gl. (10) 
Insbesondere folgt aus +1 daß ist. D. h.: Liegen die charakteristischen 

on Zahlen « von DA in der komplexen Zahlenebene im Kreise mit dem Mittelpunkte 1 und 

die dem Radius 1, so genügen die eharakteristischen Zahlen von (E-+- der Ungleichung 


u=/—1 oder /=u-l. 


3. Wir fassen die Resultate des Absatzes 2 in folgenden Satz’) zusammen: 
Satz I. Die Vektoren des Iterationsverfahrens 
Dr’); gü+)=(E+DA)ı® — Dr 


unit nichtverschwindender Determinante von D konvergieren gegen die Lösung x des Glei- 
ehungssystems Ax mit niehtverschwindender Determinante, 


I. wenn | ist, d. h. wenn die absolut größte charakteristische Zahl der 


Matrix E-+-D 4 absolut genommen kleiner als eins ist, 
(il). 
od er 


2. wenn die charakteristischen Zahlen von DA in der komplexen Zahlenebene ganz 
da) in dem Kreise mit dem Mittelpunkt 1 und dem Radius I liegen. 

Bemerkung: Ist D gleich der negativen Reziproken von A, so ıst DA=  E, und 
 charakteristischen Zahlen von D A sind sämtlich gleich — 1. Ist D also der negativen Rezi- 
ken von A hinreichend benachbart, so werden aus Stetigkeitsgründen die charakte- 
schen Zahlen von D A in der komplexen Zahlenebene noch ganz in dem Kreise mit dem 
«us I um den Punkt —1 liegen. Wir können daher dem zweiten Teile des Satzes | 
„ende anschauliche Regel entnehmen: 


Zusatz zu Satz I: Das Iterationsverfahren des Satzes I konvergiert gegen die gesuchte 

Deter- ng, wenn die Matrix D der negativen Reziproken von A „hinreichend benachbart* ist. 
Dieser Zusatz wird uns im folgenden ein Wegweiser bei der Suche nach geeigneten 
'nızen D für das Iterationsverfahren des Satzes I sein, und mit dem ersten Teil des Satzes | 
von en wir untersuchen, ob das Iterationsverfahren mit diesen Matrizen wirklich konvergiert. 


Der Satz 5 auf Seite 68 der in der Einleitung zitierten Arbeit von R. v. Mises und H. Pollnezek-Gei- 
" läßt sich als Spezialfall des Satzes I für D = const E auffassen. 


_ 
_ 
für 
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$ 2. Hinreichende Konvergenzbedingungen für das gewöhnliche Iterationsverfahren 
für lineare Gleichungssysteme. 


Hat die Matrix A eine „überwiegende Hauptdiagonale“, so ist zu vermuten, ©. 
die Reziproke der Diagonalmatrix 


A 


11 


der Reziproken der Matrix | im Sinne des Zusatzes zu Satz I hinreichend benachbart st, 
daß also das Iterationsverfahren des Satzes I mit 


D 


) 


konvergeiert. Dies Iterationsverfahren ist gleich dem Verfahren (0) von R. v. Mises uni 


H. Pollaezek-Geiringer mit €; 2 das wir, da es das gebräuchlichste ist, als «as 
„zewöhnliche Iterationsverfahren* bezeichnen wollen. Die Untersuchung ergibt: 


Satz Il.1. Wenn in der Matrix A die Diagonalglieder derart überwiegen, daß 


| v \ 
1) Max | . \ )< 1 
ik 
k=]1 
oder 
I») | \ \<1 
j 
oder 
7 
\'/a;r\ 
i,k l 


ist, so konvergiert das Iterationsverfahren des Satzes I mit 


0 
D=t 
| 
0 Ayın 
Beweis: Es ıst 
0.- 
Ayı 0 
nn nn 
Für ergibt sich daher nach (S) die Abschätzung 
n 11, 
\2 
;E+CA| „<| \ | 
MAN 
— 
i,k 1 


| | 
| 0 

( 

izek 


chr. f.aı 
th. und N 


Iren 


nuten, 


ehbart 


ISes und 


als (das 


aß 
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‚un zufolge e) in Satz II,1 die rechte Seite kleiner als I, so ist auch die linke Seite 
ser als 1 und das fragliche Iterationsverfahren konvergiert auf Grund des Satzes 1,1. 
‚\eiceher Weise ergibt sich die Konvergenz des Iterationsverfahrens aus der Ungleichung a) 
- bJ, wenn man statt (S) die Abschätzung (6) oder (7) benutzt. 
Die Ungleichungen b) und e) des Satzes II, 1 sind im wesentlichen die beiden in der 
nes zitierten Arbeit von R. v. Mises und H. Pollaezek-Geiringer (auf anderem 
\ ze) gefundenen Konvergenzbedingungen. Dort ist b) zunächst in einer etwas allgemeineren 
\orın ausgesprochen, deren Beweis jedoch auch nach unserer Methode keine neue Schwierig- 
ion gemacht hätte. 


N 


$ 3. Fehlerabschätzungen für das gewöhnliche Iterationsverfahren. 
Satz Il,2. Bei dem Iterationsverfahren des Satzes I mit 


l 0 


0 


nn 


oilt folgende Fehlerabschätzung für die i-te Näherungslösung x in bezug auf die wahre 
Lösung X: 


— 
k—1,...,n l—u, 
Ilierin ist 
n 
Max | \ ih ki Max \ ih ki ) 
| k 
und 
/ 
i,k=]1 


Beweis: Aus (11) und (12) folgt 


— -(E+ DAJi+!Dr- (E+DA)i+?Dr— , 
(E+DAy+!Dr, 
—= — 
Sun ergibt sich aus (11), daß 
(E+-DA)Jit!Dr. 


> Setzen wir dies in (15) ein, so erhalten wir 


— WW). 22 AAN. 


| Ihilfe der Ungleichung (1) folgt hieraus 


"min 
\Wurzel im Nenner können wir bei Beachtung von D=- € mit Hilfe der Ungleichung (3) 
" (4) weiter abschätzen, worauf sich die Behauptung leicht ergibt. 


*) Die Formel (15), die hier unter der Voraussetzung abgeleitet ist, daß die Reihe (12) konvergiert, gilt auch 
diese Voraussetzung. Man verifiziert dies, indem man für x(!) den Wert aus (11) einsetzt, den man noch nach 
Summenforme! für eine endliche geometrische Reihe summiert. 


20 
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Satz 11,3. Mit den Bezeichnungen des Satzes 11,2 ergibt sich die Abschätzung: 


k=1,...,n l- ur 
Beweis: Aus (16) ergibt sich 
rd) — -(E+DA4)iDr. 
Setzen wir dies in (15) ein, so erhalten wir | 
— [DAJ'(E+- DA), — . 2 2 222.019 


Durch zweimalige Anwendung der Ungleichung (1) folgt hieraus: 


+ (E+ 2 
(04 


"min 


Y yü 1) 


Hierin beachten wir. daß D— € ist und schätzen den Zähler des Bruches mit Hilfe von (5\ 
bzw. (S) ab und den Nenner mit (4) bzw. (3). Dadurch erhalten wir die Ungleichung des 
Satzes II,3 mit r=1 bzw. r=2. 


s4. Vergleich der Abschätzungen des $ 3 mit einer früheren und Durchrechnung eines Beispiels. 


I. Die Abschätzungen des Satzes II, 2 unterscheiden sich prinzipiell von der Abschätzung 
n 


in der oben zitierten Arbeit dadurch, daß bei uns die Größe Ir —ı) = (I (ap we, d": 
n 
abgeschätzt wird und dort der Ausdruck — |. 
k—1 


Wir wollen die Abschätzungen unseres Satzes Il,2 mit der Abschätzung der zitierten 
Arbeit dadurch vergleichen, daß wir mit allen Dreien eine obere Schranke der absoluten 
Beträge der Komponenten des Fehlervektors (x — xD) berechnen. Die Abschätzung der zitierten 
Arbeit ergibt 


[7] 
f v \ 
IN l IR | 
” 


Als Beispiel wählen wir das Gleiehunessvstem auf S. 65 der zitierten Arbeit (einige De 
| g 
zeichnungen sind geändert): 


32, +0,12, —0,092,=6, 0082, +4,—-0162,=12, 0052, —032, +58, = 20 


mit folgenden Näherungen des gewöhnlichen Iterationsverfahrens: 


0. Näherung = x”: 2 3 4 
1. Näherung = 1.97 3,12 4,16 
2, Näherung = x": 1,0655 3,127 4,1675 
3. Näherung = 1.096868 3,12732 4,16793 
+. Näherung = x: 1,96867 3,1274 4,1675 


2. Wir bereehnen zunächst die « und schreiben uns dazu die „durehdividierte* Ma\ 


| ) auf: 


0.05 0.03 
0.02 0,04 
001 0.06 | 


und 


von (») 
unge «les 


3eispiels. 


-hätzung 


(3)\2\1;, 


zitierten 
bsoluten 
zitierten 


(20) 


‚(20a 


De 


—Y) 


“ Matrix 
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a) Berechnung von u nach Formel (20a): 


n 
Die Zahlen \ Tih ‚“=1,2,3), sind 0,08; 0,06: 0,07. «ist ihr Maximum, also ist 
mer 


VOUS. 


b) Berechnung von «u, nach Formel (13)°): 


n 
Die Zahlen > =1,2,3), sind 0,09; 0,17: 0,12. Ihr Maximum ist 0,17. 
di 

n 
Die Zahlen 1,2,3), lauten 0,07: 0,05: 0,06. Ihr Maximum ist 0,07. 

i=-Ek 


0,17 +0,07 


Dann ist u, = —0,12. 
e) Berechnung von «, nach Formel (14): 


u,—= +3? +4? +2? +1? + 6°)'2.10-° = 0,096. 


3. Wir schätzen den Fehler der 3. Näherung ab: 


a) Benutzung der Abschätzung der zitierten Arbeit (siehe Formel (20)): 


(1+2-+2)-10°° 


55:10 


M AN = 


bh) Benutzung der ersten Abschätzung des Satzes 11,2: 


Max | 012 < 35 10° 5 


e) Benutzung der zweiten Abschätzung des Satzes 11,2: 


Max — | 10.096 <34:.10"°. 


I. kin Nachteil aller bisher ‘angeführten Abschätzungen ist, daß man mit ihrer Hilfe 
nach Berechnung der 4. Näherung nur den Fehler der 3. (nicht der 4.) abschätzen kann. 
Der Satz II,3 ermöglicht dagegen unter den gleichen Bedingungen auch noch die Abschätzung 
des Fehlers der 4. (letzten) Näherung (wir setzen r=2 in der Formel (18)): 


0,096 


10, 


S 5. 


S Konvergenzbedingungen und Fehlerabschätzungen für das Iterationsverfahren von Hertwig. 


I. Es liegt nahe, den Gedankengang am Anfang des 2. Paragraphen in folgender Weise 
zubauen: Im 2, Paragraphen war als Matrix D im Iterationsverfahren des Satzes I die 
Reziproke € der Diagonal-Matrix 


0 


>) Bei größeren Matrizen ist es zur Bereehnung von #, praktisch, sieh die beiden symmetrischen Matrizen 
dj]; Aa); 
und ik _ 


aufzuschreiben. Vgl. Dissertation, 8. 41. 
ak k d; aß IR 
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eewählt. Eine der Matrix - A"' im allgemeinen noch „besser benachbarte* Matrix erhie 


man 2. B. für » 6 in der negativen Reziproken von 


0.00 
9 


0 000 aa 


22 8 | a, 


65 ‘6 


die auch leieht zu berechnen ist. Mit ihr als Matrix D würde das Iterationsverfahren «des 
Satzes I im allgemeinen wahrscheinlich noch besser konvergieren als mit DC. 


2. Ein Verfahren solcher Art ist nun tatsächlich schon angegeben worden, und zwar von 
Hertwig”. Um zunächst zu zeigen, daß das Hertwigsche Verfahren sich als Spezialfall 
unseres Verfahrens in $ 1 auffassen läßt, übersetzen wir es in die Matrizenschreibweise: 


Gegeben sei das Gleichungssystem 


Aus der Matrix ‚I bildet man eine einfachere Matrix B, indem man in A besonders kleine 
Koeffizienten wegläßt, derart, daß das obige Gleichungssystem in mehrere für sieh eindeutig 
lösbare Gleiehungssvsteme zerfällt. Man will nun die Lösung x von (21) in Form einer 
unendlichen Reihe darstellen: 


+, +%+::.+... 


Die x; werden wie folgt bestimmt: 


Wir nehmen nun einige Umformungen vor: 
Setzen wir noch 


so wird 

=-[E+(E— A) (E— +... +(E— B’A)il 
und wir sehen aus (11), daß das Hertwigsehe Iterationsverfahren aus unserem allgemeinen 
Verfahren des Satzes I hervorgeht, wenn wir dort 


setzen. Nach Satz I konvergiert also das Hertwigsche Iterationsverfahren, wenn 


I3E—-B-14 <1 


ist. Das Erfülltsein dieser Konvergenzbedingung können wir 2. B. mit den Ungleichungen (5). 
(7) oder (S) numerisch nachprüfen. 


3. Berücksichtigen wir in (17) die Beziehung (22) und setzen 


so erhalten wir 
k=i+l 


Festsehrift Heinrich Müller-Breslan 1912, S. 37 bis 5%. 


»«) Die unter dieser Voraussetzung abgeleiteten Sätze III lauten bei anderen Anfangsbedingungen x" genaı 
vgl. aueh Anm. 


»b, Der Buchstabe C hat also von jetzt ab nieht mehr die Bedeutung, die er in den I ? und 3 hatte. 


_ 
I, = 
AL), 
| 
| 
| 
Ä 


-f. ang 
und M: 


— 


erhie 


ren «os 


war von 
ezialfal] 


(21). 
s kleine 


indeutige 
n einer 


»)\ 


meinen 


tv 
- 
.. 


(b), 


genaı 
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‚ch zweimalige Anwendung der Ungleichung (1) erhalten wir hieraus 
(An ) 
4. In gleicher Weise erhalten wir aus Gl. (19), wenn wir noch berücksichtigen, daß zu- 
(23) und (24) 


(26). 


„P'B 


5. Ganz ähnlich wie in Absatz 4 können wir auch folgende Ungleiehung erhalten 


|- 


i 
"min, 
Ilieraus ergibt sich, daß das Hertwigsche Verfahren konvergiert, wenn 


‚R’B 


"max "min 
ist. 
6. Wir fassen die Ergebnisse der ersten 5 Absätze in folgende 2 Sätze zusammen: 
Satz IIl,1. Das Hertwigsche Iterationsverfahren 
- +. tu]: 
rbr-.- 
mit DB 4-C und Det. B==0 konvergiert gegen die Lösung des Gleichungssystems Axor 


it Det. 4==-0, wenn eine der folgenden Ungleichungen besteht: 


oder 


max mım 


Satz 111,2. Für das Hertwigsche Iterationsverfahren mit den Bezeichnungen wie 
in Satz III, 1 und der Abkürzung 
2 


selten folgende Abschätzungen: 


1/, 
"min "min 
"min 


Zusatz zu den Sätzen III, 1 und 111,2: 


Sind A und B reelle symmetrische Matrizen, so kann man in den Ungleichungen 
erischen Abschätzungen wird man in den obigen Formeln die charakteristischen Zahlen 
tech ihre in der Einleitung aufgeführten Abschätzungen ersetzen. 


usw, Bei 


Imin» Max 


Wir stellen noch einen wichtigen Spezialfall heraus: 


Satz IIl,3. Das in Satz III, 1 angeführte Hertwigsche Iterationsverfahren konver- 
't für reelle symmetrische Matrizen A= (a;,) und (b;;) mit niehtversehwindenden 
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Determinanten gegen die gesuchte Lösung x, wenn y<P ist. Es gilt dann für die öte Näheru: 
l; 


Hierin ist gesetzt 


C=A-B, 
i 


$S 6. Weitere mögliche Iterationsverfahren. 


I. In manehen praktischen Fällen bietet sieh noch auf eine andere als die bisher behandelte 
Weise eine Matrix, die der negativen Reziproken der Koeffizientenmatrix des gegebenen 
Gleiehungssystems im Sinne des Zusatzes zu Satz I „hinreichend benachbart* ist. Dies« 
kann dann als Matrix D in dem Iterationsverfahren des Satzes I verwandt werden: In der 
Baustatik liegt nämlich gelegentlich zu häufig vorkommenden Belastungsaufgaben (z. B. Lösuns« 
von linearen Gleichungssystemen bei der Berechnung des durchlaufenden Trägers auf mehreren 
Stützen) zu einer gegebenen Matrix die inverse Matrix fertig berechnet vor, damit man mit 
ihrer IHilfe bei beliebig vorgegebener rechter Seite des zu lösenden Gleichungssystems_ die 
Größe der Unbekannten rasch berechnen kann. Ist nun die Genauigkeit dieses „Zahlenrecht- 
eeks" (wie es dort genannt wird) nicht ausreichend, oder will man ein System berechnen, 
das dem vorliegenden berechneten benachbart ist, so kann man zu dieser Berechnung das 
Iterationsverfahren des Satzes I benutzen. Für D ist dabei das Negative des vorhandenen 
„Zahlenrechteeks* zu nehmen °*°). 


>. Eine ganz andersartige, geeignete Matrix D finden wir durch folgende Überlegung: 
Die Matrix 4° .4 hat bei niehtverschwindender Det. A nur positive charakteristische Zahlen. 
Die charakteristischen Zahlen der Matrix e- A’ A mit e<O liegen daher in der komplexen 
/Zahlenebene sämtlich auf dem negativen Teile der reellen Achse. Wird nun außerdem 


°\<,,p., gewählt, so liegen alle charakteristischen Zahlen von e- A’ A in der komplexen Zahlen- 
"max 

ebene ganz in dem Kreise mit dem Radius Eins und dem Mittelpunkte — 1. Zur geeigneten 


Wahl von e können wir die bekannten Abschätzungen von ofen nach oben benutzen (siehe 


Kinleitung). Nach Satz I muß dann das dortige Iterationsverfahren mit D=e- A’ konvergieren, 
Ahnliehe Überlegungen lassen sich anstellen für den Fall, daß die Matrix A reell 

ıIMAaX 


Somit können wir bei Beachtung der Ungleichungen (5), (6) und (8) die Sätze aussprechen: 


symmetrisch und positiv definit ist. Wir setzen dann D=e-E mit e<O und |e|< 


Satz IV,1. Das Iterationsverfahren des Satzes 1 konvergiert für D=e-A’ mit e<V 


S 
une 
k=1 
2 
oder |e|l< 
[Max ) + 
i i,k l 
Satz IV, 2. Das Iterationsverfahren des Satzes I konvergiert bei allen reellen sy 
») 
metrischen, definiten'®) Matrizen A für D=c-E mit \e'< n 
Max 


k=] 


und e<0 für positiv definite Matrizen, und e>0 für negativ definite Matrizen. 


»°) Bei der praktischen Durehführung des Iterationsverfahrens wird man in dem hier behandelten Falle a 
bosten nach dem Hertwigschen Muster (s. Satz III, 1) vorgehen, wobei man dann für B-1 das „Zahlenrechtee: 
einzusetzen hat. 

10), PD. h. sol»he Matrizen, die zu einer definiten quadratischen Form gehören. 


n 
k 
ki kei 


f. angı 
nd M: 


— 


‚heru 


andelte 


ebenen 


Diese 
In der 


Lösuns 
»hreren 
an mit 
ms die 
nreceht- 
schnen, 
ng das 
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ber die Randwertaufgabe der gleichmäßig belasteten 


Kreisplatte. 
Von .J. Barta in Budapest. 

Inhaltsangabe: Es wird bewiesen, daß das im Titel bezeichnete zweidimensionale Problem 
hei beliebiger Einspannung auf das analoge eindimensionale Problem, d. h. auf die 
\wertaufgabe des Balkens zurückführen läßt’). Die Formel (3) besagt, daß die Durch- 
nz der gleichmäßig belasteten Kreisplatte gleich dem Mittelwert aus den Durchbiegungen 
Schar von Balken 

Bezeichnungen: 


© Durchbiegung der gleiehmäßig belasteten Kreisplatte, 
7» Durchbiegung des gleichmäßig belasteten unendlich breiten Balkens, 


p Belastung der Kreisplatte pro Flächeneinheit, 


N - p dieselbe für den unendlieh breiten Balken, 
Radius der Kreisplatte, 
d— konst. Dieke der Kreisplatte oder die Dieke des unendlich breiten Balkens, 
 Blastizitätsmodul, 
(uerkontraktionszahl, 
Eo' 


Plattenkonstante. 


Die Differentialgleichung für die Durchbiegung der dünnen Platte lautet 


Wir wollen uns mit derjenigen Lösung [ befassen, die am Kreisumfange die Randwerte 


(erste Randbedingung) 
und 


7 (zweite Randbedingung) 
| 


annimmt. #2 ist die nach dem Innern des Kreises gerichtete Normale. Z und Z sind gegebene 
Kunktionen des Kreisbogens. 

Irgendein Punkt der Kreisfläche sei mit bezeichnet (Abb. I). Die Endpunkte der 


Sehne dureh den Punkt P seien mit 0, 0' und die entsprechenden Randwerte mit [, < 
bezeichnet. Die Normalen » und bezeichnen die Riehtungen von @ nach bzw. 


AS93S1 


Abb. 1. Abb. >. 


') Zur Lösung der Plattenprobleme, welehe in der Baustatik vorkommen, benutzen die Ingenieure oft solehe 
den, die das Plattenproblem auf das Balkenproblem zurückführen. Aber diese Methoden liefern nur mehr oder 
srobe Annäherung der genauen Lösung des Problems. Man findet eine knappe Übersieht über diese Methoden 
ı: Gehler, W.: Balkenbrücken. (Handbuch für Eisenbetonbau.) 3. Aufl., S. 320. 


| N 

| 
| 
= \ | | 
Le | 

| 
‚chtecl | 
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von 9' nach ©. Die Riehtungen von Q nach P bzw. von Q’' nach P seien mit m und 


bezeichnet. Der Wert : = für den Punkt © und yon, für den Punkt 9 ist durch die Ra. 


bedineungen bestimmt. 
Stellen wir uns einen Balken zwischen Q und 0 vor. Der Balken habe die Stützw« t« 


= Er sei unendlich breit und von konstanter Dieke. Seine Belastı 


- - 


- 


habe «den bezogenen Wert 4=75 pP. Die Differentialgleichung für die Durehbiegung dieses 
Balkens lautet 
eng 


Als Randbedingungen für den Balken seien die Randbedingungen der Kreisplatte vor: 
veschrieben. Es sei also 


dm om 


für den Punkt © und 


dy 0% 


dm dm 
für den Punkt @. 


Jeder Sehne, welehe dureh den Punkt P°_ geht (Abb. 2), entspricht nach unsrer Vor- 
schrift ein Balken, dessen Durehbiegung im Punkte P wir mit np bezeichnen. Auf diese 
Weise haben wir uns Balken von unendlich großer Anzahl vorgestellt. Wir bezeichnen den 
kleinen Winkel zwischen benachbarten Sehnen mit dgy. Die Durcehbiegung der Kreisplatte 
im Punkte P läßt sich dann durch die Formel 


IT 


0 


ausdrücken. Unser Plattenproblem ist also durch diese Formel auf das Balkenproblem zurück: 
oeführt. Wir wollen im folgenden die Richtigkeit der Formel (3) beweisen. 


/um Beweis schreiben wir die Lösung der Gl. (1) in der Gestalt 


d 
7 
»2 29 | | R:? 19 Q n 
0 
auf: das erste Glied 
ist die Lösung der Randwertaufgabe I1.1Z, : 7 für den Kreis mit den Randbedingungen 
| 


Das zweite Glied der rechten Seite der Gl. (4) ist die Lösung der biharmonischen Randwt 
aufgabe 1.1/, 0 für den Kreis’) mit den Randbedingungen 


N 
_ Go 
ud - — 


on dm’ 


Sieh» z.B. das Bueh: Riemann-Weber: Die partiellen Differentialgleichungen der mathematischen Phy 
Ss, Aufl. von Frank und Mises. 1. Teil, S. 851, Formel (21). 
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Ist s— Ra die Bogenlänge, so haben wir 


dn 
om oOndm Osdm’ 
hei 


OL 


An 
dm 


Es ist mithin 


| 
_to(V0OP- 
Wir haben aus dem Dreieck den Winkel 


nach dem Sinussatze 


I 


Diese Gleichung drückt den Zusammenhang zwischen und 
Durch Differentiation der letzten Gleichung erhalten wir 


sin(p  sin(p a) 


1 
cos(y cos (q da). 
bzw. 
Rcos(y a) 


Aus dem Dreieck 0 PO der Abb. 1 haben wir 
o=yR?+r?—2Rrcosia- 
Die bekannte Kreiseigenschaft o 0’ = — ergibt 
PP? 


oder 
ar? 
- R?—r 


TyR+r—2Rrcosia— 
Aus der Abb. 1 haben wir noch 

cos 


Mit Rücksicht auf (7) und (8) 
R- reosla 


IT — 


folgt noch 
rsin(ia 


rcosia 


tg (09 R 


Die Gl. (5) sehreibt sich also 


05 VYR+r—2R recos(a—d) , dc rsın la 
on om rcosia — 
ragen wir diesen Ausdruck in (4) ein, so ergibt sich 
r sin (a — | 2I[R 


und 


ds 
> — sın - 


(D). 


aus, wenn der Punkt P fest ist. 


(6). 


(S). 
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Wir formen das’ zweite Integral in der letzten Gleichung durch die partielle Integration 


de) 
da 
da 


M 


um, wobei wir der Kürze wegen die Bezeichnung 
rsin(a 


9 - 
R|R r cos (a 


Nach Durchführung dieser Rechnung ergibt sich 


benutzen. 
en 
Ä 
-2Rrcosla— 
- 6) > da 
Aus der Gl. (2) folgt dureh Integration 
+Aw+BwW+Cu+D (10) 
und 
dn 
du 6X + 
Wenn der positive Ast der Koordinatenachse « mit PQ zusammenfällt und zu «= 0 der 
Punkt 7’ gehört, so ist aus (10) 
p=D (il). 
Aus den Randbedingungen 
ergeben sich die Gleiehungen 
34 34 N? Ao’+Beo’-Co+D=/, 
— 3A0o?—-2Bo—Ü = = 
Dureh Auflösung dieser vier Gleichungen lassen sich die Konstanten A, B,C', D berechnen. 
So erhalten wir 
Nach (11) können wir schreiben 
de \ =— ” —— + 


7 
"Id 00° -300?+ 0° 


oder 
d 


Die rechte Seite dieser Gleichung geht, wenn man die Formeln (6), (7), ($) und q= 8 pberüc 
sichtigt, in die rechte Seite der Gl. (9 über, wodurch die Richtigkeit der Formel (3) © 


wiesen Ist. 


.f.angı 
und Me: 


tion 


la 


(9). 


(10) 


- () der 


(il). 


m 


ecehnen. 


berüc! 


(3) © 
59 
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KLEINE MITTEILUNGEN 


„um Newtonschen Näherungsverfahren. 
Newtonschen Näherungsverfahren geht 
hekanntlieh von einem Näherungswert aus, für 
Funktionswert und zweite Ableitung das 
"e Vorzeichen haben. Man erhält dann den 
‚oJwert mit genürender Genauigkeit durch um 
enieer Schritte, je geringer die Krümmung 
;wischen Näherungswert und  Wurzelwert 
nden Kurvenstückes ist. Es liegt daher nahe, 
ııf der X-Achse statt der gleichmäßigen Skala 
"unktionsskala aufgetragen zu denken, die 
hetreffende Bogenstück in passender Weise 
\t, Dabei muß man aber zweierlei beachten: 
ı\ darf die Streekung nieht so stark sein, daß 
Krümmung der Kurve ihr Vorzeichen ändert, 
‚weitens muß die die neue Skala bestimmende 

‚ktion so einfach sein, daß die Arbeit für die 

rische Auswertunz der neuen Näherungs- 
») nieht wesentlich größer ist als die bei der 

söhnliehen Newtonschen Formel. 

>». Es seien die Nullstellen der Funktion 


bestimmen: den n-ten Näherungswert be- 
rechnet man gewöhnlich aus dem (n — 1)-ten 7,__ı 
der Formel 

r r (N 

wo (dureh den Strieh auch weiterhin die Ableitung 
nach .r bezeichnet werden soll. Denken wir uns 
nun auf der X-Achse eine Funktionsskala 


aufeetraeen und als Ordinate dazu unsere obige 


Funktion 


un auf diese Kurve von dem gleichen Näherungs- 
wert (f„_ı) aus das Newtonsche Ver- 
Iualıren angewendet, so erhält man einen anderen 
n-ten Näherungswert r, aus 


(zu-ı) - 
dz 


dr 
d z 
ir die Reehnung wird in dieser Formel wieder z 
tureh 7 ersetzt. Die Formel (1) läßt sich nun auch 
‘olzendermaßen schreiben: 


r r 
un = 
) 
1 
In-ı) 
n nA 


Iep 


us sieht man folgendes: Ist der Differenzen- 
zwischen .r„, und für die transfor- 
ende Funktion kleiner als ihre Ableitung für 
den (n —1)-ten Näherungswert und haben, was im 
| einen zutreffen wird, beide das gleiche Vor- 
eichen, ist 


— In-ı| > 
ie Anwendung der Newtonschen Formel 


° transformierte Funktion führt unter der 
‘emachten Voraussetzung schneller zum Ziel 


als die Anwendung auf die ursprüngliche Funk- 
tion. Haben .r und 9’ (r) das zleiche Vorzeichen 
und ist 9” (r) in dem in Frage kommenden Gebiet 
positiv, ist also die g (r) darstellende Kurve nach 
oben konkav, so ist bei einer auf den Nullpunkt 
zu laufenden Näherungsfolge der Differenzen- 
auotient kleiner als die Ableitung im — 1)-ten 
Näherungswert, die Konverzenz bei der trans- 
formierten Funktion ist daher besser als bei der 
Ausgangsfunktion. Bei nach oben konvexer Funk- 
tion $ (r) wird dazeren in diesem Fall die Kon- 
vergenz der sich vom Nullpunkt entfernenden 
Näherungsfoleen durch Transformation ver- 
hessert. Haben und (r) entzerenzesetztes 
Vorzeichen, ist es umzekehrt. Geht man also 
vonzu großen Näherungswerten aus, 
kann man dureh Transformation mit 
"Funktionen, deren zweite Ableitung 
indeminBetrachtkommenden Gebiet 
positivbzw. negativist.jenachdem, 
ob.rundg(r) zleichesoderentgeren- 
vesetztes Vorzeichenhaben.dieKon- 
vegen von zu kleinen Näherungs- 
werten aus, muß die zweite Ablei- 
tunze der transformierenden Funk- 
tion das entzerenzesetzte Vorzei- 
chenhaben,. wieinobiren Fällen. 
bei ist aber darauf zu achten, daß die zweite Ab- 
leitung der ursprünglichen und der transformierten 
Funktion das gleiche Vorzeichen. nämlich das der 
Funktion an dieser Stelle haben. 


3. Reehner'sch am einfachsten werden im allre- 
meinen die Formeln, wenn man (ir) x" setzt, 
wobei wir uns auf positive x beschränken. und 
zwar muß. da hier z und g’ (r) das gleiche Vor- 
zeichen haben, r>1 sein. wenn man von einem 
zu großen. r<1. wenn man von einem zu 
kleinen Näherungswert auszeht. Man erhält dann 
die Näherungsforme! 


r 
In >= / In 
Für eine trinomische Gleichung 2" + + B= 
erhält man z. B. 


(D). 


/ (m —r) Am-ıP —rB 
27 | 
dabei müssen die in Betracht kommenden r in 
einem Gebiet liegen, in dem 


bzw. 

m(m—r) em -?2r LAp(p— r)am—!r (7) 
das zleiche Vorzeichen wie die Funktion haben. 
Die Formeln (5) und (7) werden besonders einfach. 
wenn mn=-r oder p=*r ist, da dann der Zähler 
von (5) aus nur zwei Summanden besteht. und (7) 
nur ein Glied hat. Nehmen wir etwa die Gleichung 


0,162 


die die Wurzeln — 0,24: + 0,4; + 0,6 hat, so kann 
man als Näherungswert für die größte Wurzel 
x 1 nehmen, einen Wert, für den Ordinate und 
zweite Ableitung positiv sind. (7) wird in dem in 
Frage kommenden Gebiet für r=2 und r=3 
ebenfalls positiv. Für r 1 erhält man aus der 
Newtonschen Näherungsformel 


2 + A2u-,°--B 
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die Näherungsfolze 1: 0.799: 0,681: 0,623 . . . 
r=?2 ergibt sieh aus 


die Folze 1: 0,773: 0.659; 


aus 
En 


die Folge 1: 0,734: 0,633: 0,602 . . 
etwa die Gleichung 


Für 


—2B (10) 


0,611... und fürr=3 


3 
A 


B2._ı +24 


(11) 
Wählen wir 


— 38.09 x? —- 3.3124 =0 (12) 
mit den drei reellen Wurzeln — 0,91; + 1.4: + 2.6. 
so kann man als Näherungswert für die mittlere 
Wurzelr -1.2wählen. Funktionswert undzweiteÄb- 
leitung haben hier das gleiche positive Vorzeichen. 
7) ist selbstverständlich für r = 1, aber auch für 


1 
r positiv, dagegen für r 3 schon 


Aus der Formel (9) berechnet man mit =1,2 


negativ. 


die nächste Näherung rs = 1,389; für r= _ erhält 


man aus 


+ Am-ı?- (13) 


| 
n | 
n—i | 


mit 


oeibt für r 


| 

mit 2,=12 7, =1,401. r, liest hier also schon 
auf der anderen Seite des Wurzelwertes, wie das 
bei Ausgang von einem Kurvenpunkt sein kann, 
in dem Funktionswert und zweite Ableitung nicht 
mehr das gleiche Zeichen haben. 


4. Benutzt man für die Transformation die Funk- 
tionen 


2 bzw. 2 In.r ?), 


so erhält man die Formeln 
und 
(xn-ı) 
. 


die erste konvergiert besser bei auf den Nullpunkt 
zu, die zweite bei vom Nullpunkt fortlaufenden 
Folgen. Für die Gl. (8) haben Funktionswert und 
zweite Ableitung der transformierten Funktion in 
dem in Betracht kommenden Gebiet das zleiche 
Vorzeiehen: Formel (15) darf also angewandt wer- 
den und liefert mit .r, ==1 als nächste Näherung 
2, = 0,776, also besser als die Newtonsche aber 
in diesem Falle nieht so gut wie Formein (10) und 
(11). Für Gl. (12) haben Funktionswert und zweite 
Ableitung der transformierten Formel für .r, = 1.2 
entzexengesetztes Zeichen. Wendet man trotzdem 
Formel (16) an. erhält man mit 7, = 1.2 den Wert 
1,45, also einen Wert jenseits des Wurzel- 
wertes. 


Freiberz (Sa.). Fr. A. Willers. 612 
I!) Beier: Unterriehtsblätter 4? (1936). S. 9. 
Bailey: Phil Mag. (7) 18 (10%4), 8.529 bis 539 setzt 
°—loe x, um bei der Berechnung der Nullstellen ganzer rat. 
Funktionen bequem Logarithmen anwenden zu können. 


Math. und 


Minimaisinkgeschwindigkeit von 
zeugen. \Wenn 


b die Spannweite, 
F die Flügeltläche, 
.„—b?:F das Seitenverhältnis, 
c ; den Auftriebsbeiwert des Tragflügels. 
Cıwp den Beiwert des Profilwiderstandes (les 
Tragflügels, 
Cırsj den Widerstandsbeiwert eines nichttraxn. 
den Flugzeugteiles, 
Fj die Fläche des betreffenden nichttragen. 
den Flugzeugteiles, 
z den Ausdruck I Fj summiert über 
alle niehttragenden Flugzeusteile, 
( das Fluggewicht eines Flugzeuges, 
die Beschleunigung der Erdschwere, 
y das spezifische Gewicht der Luft 


bedeutet, dann gilt für die Sinkgeschwindigkeit 
eines Flugzeuges 


( 

29. 

wenn elliptische Auftriebsverteilung vorausgesetzt 
wird. 

Der Zusammenhang zwischen und liegt 
auf Grund der Anblasversuche derart vor, dab ec, 
als Funktion von cwp angesehen werden kann. 
Für eine vorgegebene Spannweite und konstante 
Werte von G,g.y und 2 ist dann v, eine Funktion 
von und Cwp, die einen Extremwert!) für jene 


Werte cıwrpo und /, hat, die die Gleiehungen 
Ö 
und -—0 
O4 
befriedigen. 


Die Durchführung der Differentiation liefert die 
Gleichungen 


32 


32 
Be, 3Cwpol. — ( ==), 
dc. 
wobei = . genommen an der Stelle cr 
tCwp 


bedeutet. Multipliziert man die erste Gleichung 
mit €’ po und addiert sie zur zweiten, so erhält man 
als bedeutsames Ergebnis 


„10 
Eo; 
Cw po 
wobei E, den maximalen Gleitwert (reziproket 


Gleitzahl) des verwendeten Tragtlügelschnittes Für 
— eroßes Seitenverhältnis darstellt. 

Die weitere Rechnung ergibt für das günstigste 
Seitenverhältnis 


a Ci Po 
Av 


Für sehr große Spannweiten geht dieser Ausdruck. 
nach binomischer Entwicklung der Quadratwur 
über in 

CH 

40 . 
po 


Eine Überschreitung dieses Grenzseitenverhältni> > 


ergibt keinen Gewinn hinsichtlich Sinkgeschwint 
keit. 


1) Man überzeugt sieh leicht, daß dieser Extremwer 
Minimum ist! 


1.2 len Wert = 1.398 und schließlich 


r.f. ang w., 16, Heft 5 


2917 
‚und “bar 1936 Buchbesprechungen 317 


n Flug. man mit den Werten Cwpo und 4, in Eine eingehende Diskussion der Faktoren und » 
beziehung für die Sinkgeschwindigkeit ein, im Hinbliek auf die Abhängigkeit von E, ergibt: 
't man, wenn man zur Abkürzung Für größer werdendes EZ, wird «u, somit die Sink- 
geschwindiekeit kleiner und », somit die Flugge- 
/ schwindiekeit größer. 
Fir 
lürels Zusammenfassend kann folgendes festgestellt 
ugels, - 
indes (les werden: 


Für Flugzeuge. die bei kleinster Sinkgeschwindig- 


httrasen.- 9, keit möglichst hohe Fluggeschwindigkeit erreichen 
80 Y sollen, sind Tragtlügelprofile mit möglichst hohem 


httragen- E,. d. h. Gleitwert bei unendlich großem Seiten- 


verhältnis zu wählen. Alle Profile mit gleichem 
E, sind hinsichtlich kleinster Sinkgeschwindigkeit 
eleichwertigz.  KEntscheidend für die Wahl eines 


liert über 


le, 


ob? 


res, 
were, 


ft 


vindiekeit 


rusgesetzt 
Cwp liegt 
r, dab ec, 
len kann. 
konstante 


Funktion 
) für jene 


ijefert die 


telle 


Gleichung 
rhält man 


(reziproke 
nittes lür 


Fünstigste 


Ausdrı K. 
ratwure 


hältni> 
chwind 


enwer: 


\nalor ergibt sich für die Fluggeschwindigkeit 


I)ie Koeffizienten @ und » hängen. was die Profil- 


oioenschaften anbelangt, nur von E,= ab, 
Cw po 

s ereeben somit unter sonst gleichen Umständen 

alle Flügelprofile mit gleichem EZ, gleiche kleinste 

Sinkzesehwindiekeit. gleiche zugehörige Flugge- 

schwindiekeit und somit gleiche zugehörige Gleit- 

zahl. 


Protiles werden die Eigenschaften bei kleinen 
Werten sein, denn auch in diesem Bereich werden 
für ein hochwertiges Flugzeug relativ kleine Sink ge- 
schwindiekeiten anzustreben sein. Das zu wählende 
eünstigste Seitenverhältnis ist bei gerebenem E, 
proportional ci po, man wird daher von allen gleich- 
wertizgen Profilen jene in engere Wahl ziehen, «ie 
auf ein Seitenverhältnis führen, bei dem das Trag- 
tlügeleewicht ein Minimum erreicht. 

Die Werte € go. po und eharakterisieren ein 
Tragtlügelprofil in weitgehendem Maße. 

Die angegebenen Formeln ermöglichen den 
direkten Entwurf von Flugzeugen kleinster Sink- 
eeschwindigekeit. 


Knittelfeld. K. Novotny. 592 


BUCHBESPRECHUNGEN 


Dr..Ine. KARL KARAS, Privat- und Honorar- 
lozent a. d. Deutsch. Techn. Hochsch. Brünn, Die 
kritischen Drehzahlen wichtiger Ro- 
ıoren. IV + 1548. m. 40 Textabb. u. 22 Zahlen 
taf, Wien 1935, Verlag Julius Springer. Preis brosch. 
I8 M. 


Die Erkenntnis von der großen technischen Be- 
leutung kritischer Zustände ist heute wohl Gemein- 
„ut aller Ingenieure geworden; aber die richtige Be- 
wehnung kritischer Drehzahlen auch dann, wenn 
keine einfachen Grundfälle vorliegen, birgt trotz 
vieler guter Lehrbücher und Formelsammlungen 
noch viele schwierige Probleme in sieh, denen der 
\erfasser des vorliegenden Buches mit der Absicht 
‚u Leibe geht. dem Konstrukteur tunlichst ohne um- 
ungreiche Studien und verläßlich die Vorausbereeh- 
nung zu ermöglichen. Ein wichtiges Hilfsmittel da- 
‚u ist das bekannte und technischen Problemen so 
zus anpassungsfähige Verfahren von Ritz. Es 
wirt in dem Buche besonders auf die kegelstumpf- 
und stufenförmigen Rotoren (von Dampfturbinen) 
inzewendet, die in der vorliegenden Schrift erst- 
nulr einer einheitlichen Bearbeitung unterworfen 
worden sind. Nebenbei ergeben sich allerlei neue 
wenige bekannte Zusammenhänge und Resul- 
10, so die Verhinderung der kritischen Drehzahlen 
eh die Kreiselwirkung der aufgekeilten Scheiben. 
'e schon von Grammel bemerkt) hier in ganz all- 
„oneiner Art in Erscheinung tritt. Für den Ver- 
“er der klassischen Mechanik ist es außerdem er- 
ch, daß auch in diesem Buch wieder einmal 
Iamiltonsche Prinzip zur Lösung tech- 
soer Fragen herangezogen wird, und daß es dabei 
universelle Bedeutung auch in diesem Spezial- 

leutlich erweisen konnte. 

s dem Inhalt sei hervorgehoben, neben der 
'sätzlichen Problemdarstellung und deren Lö- 
einschließlich von Näherungsformeln noch die 
ıllune symmetrischer Fälle von Trommeltur- 
nit kegelstumpfförmigem Läufer, weiter die 


J 


Berücksichtigung der Kreiselwirkung der Läufer- 
scheiben, Betrachtungen über das Curtis-Rad, die 
Gleichungen von Brunelli angewendet auf die all- 
gemeine Zweifelderwelle mit Kreiselwirkung und die 
Zweiflußtrommelturbine mit abzesetztem Läufer, 
endlich die Betrachtung der allgemeinen Dreifelder- 
welle ohne Kreiselwirkung. Daneben stehen noch 
einige Kapitel, die mehr allgemeinen analytischen 
Charakter besitzen, und die besonders zur schärferen 
3estimmunge der ersten und zweiten kritischen 
Drehzahlen herangezogen werden. 

Schon aus diesen kurzen Andeutungen mag der 
reiche Inhalt des Buches erkannt werden: «dazu sei 
noch auf die flüssige Darstellung sowie das reich- 
haltige Literaturverzeichnis hingewiesen. Es kann 
wohl gesagt werden, dab der wissenschaftlich arbei- 
tende Ingenieur, daneben aber auch jeder, der für 
technische Mechanik im modernen wissenschaftlichen 
Gewande Interesse besitzt und die Sprache der 
analytischen Mechanik genügend beherrscht. von 
diesem Buche reiche Anregung und Befriedigung 
empfangen wird. 


Hannover. A. Pröll. 598 


ZARISKI, A\lzebraie Surfaces. (Ergeb- 
nisse der Mathematik und ihrer Grenzgebiete, hrsg. 
v. d. Sehriftltg. d. „Zentralblatt f. Mathematik“, 
3. Bd. 5.) V1+ 198 8. Berlin 1955, Julius Springer. 
Preis 22,75 M. 


Erste und zusammenfassende und moderne Dar- 
stellung der Theorie der algebraischen Flächen, wie 
sie sich hauptsächlich in der italienischen Schule, 
der der Verfasser angehört, entwickelt hat. Dabei 
wird insbesondere auf die Theorie der mit einer 
algebraischen Fläche verbundenen transzendenten 
Funktionen eingegangen und auf die Zusammen- 
hänge mit der Topologie. 


Hamburg. 


Blaschke. 600 
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Buchbesprechungen 


Ztschr. f. anı w. 
Math. und M 


B. L. VAN DER WAERDEN, Gruppen von 
linearen Transformationen. („Ergebnisse 
der Mathematik und ihrer Grenzgebiete“, hrsg. v.d. 
Schriftleitunze des „Zentralblatt für Mathematik“, 
4. Bd. 2.) 91 8. Berlin 1935, Verlag Julius Springer, 
Preis brosch. 8,50 M. 

„Gesamtdarstellung (dieses Gebietes, die zugleich 
die Beziehungen zur Theorie der kontinuierlichen 
Gruppen und zur projektiven Geometrie klar hervor- 
treten läßt.“ 

„Die Isomorphismen (der orthogonalen Gruppen in 
densingulären Fällen » = 3, 4, 5, 6 werden hier unter 
Hervorhebung ihrer reichhaltigen geometrischen und 
algebraischen Beziehungen von Grund aus neu her- 
zeleitet.“ 

Da die Gruppentheorie für die Physik wiehtig zu 
werden beginnt, dürfte dieses treffliche Buch auch 
für manchen „angewandten“ Mathematiker brauch- 
bar sein. 


Hamburg. Blaschke. 600 


W. N. BAILEY, M. A. D. Se. Senior Leeturer in 
Mathematies in the University of Manchester, 
Generalized Hypergeometrie Series. 
Cambridge Traets in Mathematies and Mathemati- 
cal P’hysies, No. 32.) 108 8. Cambridge 1935, Uni- 
versity Press. brosch. 6/6 Sh. 

Untersuchungen über hypergeometrische Reihen, 
verallgemeinerte hypergeometrische Reihen, Hei- 
nes Reihen und Appels hypergeometrische Funk- 
tionen zweier Veränderlicher. 


Hamburg, Blaschke. 600 


Irof, Dr. WITTING, Oberstudienrat i. R. Repe- 
titorium und Aufrabensammlung zur 
Differentialreechnung. Sammlung Göschen 
146. erlin u. Leipzig 1935. Verlag Walter de 
Gruyter & Co. 139 S. mit 58 Fig. Preis 1,62 M. 

Das vorliegende Bändchen ist eine vollständige 
Neubearbeitung des früher von Junkers herausgege- 
benen  zleiehnamizen Bändcehens der Sammlung 
(Göschen, (die sich an (die vom gleichen Verfasser 
herauszerzebene Differentialreehnung anschließt. 
Kine Durcharbeitung der 405 gutgewählten Beispiele 
und Aufgaben, die den dort behandelten Stoff zum 
Teil noch erweitern, dürfte eine ausreichende Fertig- 
keit in der Handhabung der Formeln und Methoden 
der Differentialreehnunge gewähren. 


Fr. A. Willers. 577 


Freiberg (Sa.). 


Dr. ARTHUR HAAS, Prof. für Physik a. d. Uni- 
versität Wien, Atomtheorie. Vl1ll+292 8. m. 
st Fire. und 5 Tafeln. Dritte. völlie umgearb. u. 
wesentlich verm. Aufl. Berlin und Leipzig 1936, 
Verlae Walter de Gruyter u. Co. Preis geh. 8.50 M, 
veb. WM. 

Gerenüber der 192D erschienenen, in «ieser Zeit- 
schrift IN, S. 345 von v. Mises besprochenen zwei- 
ten, trägt die neue Auflage den schnellen Fortschrit- 
ten der Atomphysik durch sorgfältige Überarbeitung 
und Erweiterung Rechnung. Insbesondere sind die 


neueren Forschungen auf dem Gebiet der Kern- 
physik und der Spektroskopie und «die neueren 


Kenntnisse über Korpuskularstrahlen und über (die 
Struktur der Atomhülle berücksichtigt worden. Das 
Buch hat die bekannten Vorzüge fast aller Haas- 
schen Werke: Kurze und «doch verständliche, über- 
aus klare Darstellung auch «ler neuesten Erzebnisse, 
läßt aber wohl die noch bestehenden Schwierirkei- 
ten und Unstimmigkeiten nieht immer zenürend her- 
\vortreten. 
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Dr.-Ing. habil. F. WEINIG, Dozent an der Ten. 
Hochschule Berlin, Die Strömung um ir 
Schaufelnvon Turbomaschinen. VII -+ 
142 S. m. 120 Abb. Leipzig 1935, Verlag von J. A 
Barth. Preis geb. 17,50 M. 

Die Arbeit macht das geradlinige Schaufelgiiter 
mit beliebiger Profilform der rechnerischen Beha | 
lung dadurch zugänglich, dab die Profile dure) 
gleichwertige Streckenprofile ersetzt werden. De, 
beweis der Gleichwertigkeit letzterer wird erbracht. 
Das geradlinige Streekenprofilgitter wird zunächs: 
unter Heranziehung früherer Arbeiten des gleichen 
Verfassers durchgerechnet und das Ergebnis in 
Kurventafeln dargestellt, aus denen insbesoniler 
das Verhältnis der Auftriebsziffern eines Strecken 
protiles in Gitteranordnung zu der eines gleichen 
alleinstehenden Iroftils bequem entnommen werten 
kann. Von besonderem Wert sind die hierzu ze- 
gebenen Näherungsformeln für den Fall enger. 
Teilung (v1 <0,7), die in der Praxis recht häufig 
vorkommen. Auch die Austrittsablenkung der Strö 
mung wird ermittelt und insbesondere ihre Änderung 
in Abhängigkeit des Anströmwinkels an Hand eines 
einfachen graphischen Verfahrens gezeigt. 

Die Behandlung beliebiger Protile setzt die Er 
mittlung des gleichwertigen Streckenprofils voraus. 
d. h. seiner Länge und Neigung sowohl für (das 
Gitterprotil wie für das Einzelprotil. Dies gestalte: 
sich bei Kreisbogenprofilen einfach und ist dure) 
Kurventafeln erleichtert. Auch für beliebige Profils 
sind die Verfahren angegeben, und die Geschwinldig 
keits- und Druckverteilung längs eines Gitterproiil 
ermittelt, ebenso wie umgekehrt die Bestimmung 
eines Profils für gegebene Druck- und Geschwinldig 
keitsverteilung gezeigt ist, so daß also auch die 
kücksicht auf Vermeidung der Kavitation beachtet 
werden kann. Doch sind diese Verfahren nur bei 
dünnen, schwach gewölbten Profilen oder für Pro- 
file mäßiz weiter Teilung ohne erheblichen Zeitauf- 
wand durchführbar. Für gegebene beliebige Profile 
ist der Entwurf des Strombildes nicht zu umgehen. 

Wenn es deshalb auch fraglich ist, ob die Ergel 
nisse der vorliegenden Arbeit beim Entwurf von 
axialen Strömungsmaschinen in erheblichem Umfang 
werden Verwendung finden können, so ist es doc) 
dankbar zu begrüßen, daß der ganze schwierige Stoff 
zusammengefaßt und gezeigt ist, welcher Weg (lei 
Lösung im einzelnen Fall beschritten werden kam. 
An Beispielen ist gezeigt, dab die Ergebnisse mit 
denen anderer Verfahren sowie mit dem Versuch 
gut übereinstimmen. Anzuerkennen ist ferner (li 
große Zahl von Strömungsbildern, die das Verständ- 
nis erheblich erleichtern. Schade ist, dab anderer 
seits die Benutzung der Arbeit dadurch erschwer! 
wird, daß im Text Hinweise auf vorausgegangen« 
Ableitungen fehlen, z. B. durch Angabe frühere 
Gleichungen. Die Darstellung ist ferner trotz (ler 
sichtbaren Bemühungen des Verfassers gerade 
der wichtigen zweiten Hälfte des Buches noch \ 
sehr auf den Gedankenkreis desjenigen abgestellt. 
der sehr viel mit den Verfahren der konform 
Abbildung zu arbeiten hat. Deshalb wird sich (le 
Ingenieur, der mit dem Entwurf von Strömung» 
maschinen betraut ist, auch wenn er die Gr 
lagen dieser Verfahren kennt, vielfach schwer 
reehttinden, zudem ihn die mathematischen \) 
leitungen weniger angehen als die Ergebnisse. 
dureh dürfte die sehr wünschenswerte Verbreit "+ 
des Buches gerade in den Kreisen, für die es 0! 
bar geschrieben ist, leiden. Hoffentlich gibt 
bald Gelegenheit, diese formalen Mängel der 
verdienstvollen Arbeit zu beheben. 


C. Pfleiderer. 
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-HABERT, Wärmetechnische Tafeln 
‘erlagen für die kechnungen des Wärme- 
nieurs in >chaubildern und Zahlentafeln), 
>75 8.. 36 Tafeln. Verlag Stahleisen m. b.H., 
seldorf. u. J. Springer, Berlin 1935. Preis 14,.50M, 

ir Wasserdampf, Luft und technische Gase 
»t der Wärmeingenieur in diesem Werk eine 
mmenstellung aller erforderlichen Zahlenwerte. 

beznügt sich der Verfasser nicht allein mit 
Siehtung dieser Werte oder ihrer Umrechnung 
einheitliches Maßsystem, sondern er bringt alle 
'onkliehen Darstellungsarten,. so daß sich in 
\en Fällen auch Integrationen oder Umzeich- 
ven erübrigen. Es sind beispielsweise für 
\sserdampf 5 verschiedene Zustandsschaubilder 
jeweils anderen Koordinaten angegeben. Bei 
sen findet man neben den polytropischen Be- 
ungen auch die wirklichen Adiabaten technischer 
abgebildet. Selbstverständlich sind nicht alle 
/shlenangaben mit der gleichen Vollkommenheit 
bokannt wie die der Dampftabellen, aber gerade 
derartige Werte wie z. B. Zähigkeit von Gasen 
um! Wasserdampf oder Zustände von Gasen bei 
hohen Temperaturen sind dafür gewöhnlich nur in 
der Literatur verstreut zu finden. Das umfang- 
reiche Werk ist aus einem Bedürfnis der Praxis 
heraus entstanden und wird sich sicher einen 
latz unter den Hilfsmitteln des Wärmeinzenieurs 
erobern, der die mühevolle Arbeit des Verfassers 
lolınt. In einem Anhang ist in losem Zusammenhang 
mit dem übrigen Stoff eine Literaturübersicht von 
Dr.-Ing. H. Schwiedeben beigefügt über 
uerungstechnische Berechnungen der betrieblichen 
Praxis. Busemann. 616 


Verhandlungen der in Tallinn und Tartu vom 
20. bis 23. August 1935 abgehaltenen achten Tagung 
der Baltischen Geodätischen Kommission, redigiert 
vom Generalsekretär Imari Bonsdorff. 296 S. 
1986, Osakeyhtiö Weilin & Göös Aktiebolag. 

Diese Veröffentlichungen über die Tagungen der 
baltischen Geodätischen Kommission sind aus zwei 
(ründen wertvoll: erstens geben sie ein Bild von 
den Stand der Arbeiten in den beteiligten Ländern, 
‚weitens zeigen sie die derzeitige Entwieklungsstufe 
er in Frage kommenden Meßverfahren. Nach dem 
bericht des Generalsekretärs I. Bonsdorff neh- 
nen «lie von der Kommission beabsichtigten Arbei- 
en noch mindestens 12 Jahre in Anspruch. Die 
Grundlinienmessungen sowie die astronomischen 
ingen-, Breiten- und Azimutbestimmungen können 
‚ir «(len ganzen Ostseering als beendet angesehen 
werden. 

Neben den dynamischen Schweremessern werden 
unmer mehr jetzt auch statische Schweremesser ver- 
wendet. E. Kohlschütter teilt in einem Be- 
ehte mit, daß die Entwicklung des offenbar sehr 
statischen Schweremessers von 
I Haalck abgeschlossen ist, und daß bei diesem 
Instrument für Messungen auf dem Lande und auf 
or See verschiedene Formen nötig sind. Die Ge- 
uuirkeit der Messungen auf dem Lande betrug bei 
“ner Monatsleistung von 150 Punkten + I mgal: bei 
Messungen auf fahrenden Schiffen wurde eine Ge- 
uuliekeit von #3 bis #5 mgal erreicht. 
‚linzuweisen ist hier auf die Berichte über die 
nilinienmessungen mit 24 m langen Invar- 
en und deren Eichung sowie die Ausführung 

"“einnivellements. 

'resden, 
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ner sind bei der Schriftleitung folgende 
er eingegangen (ausführliche Besprechung 
't vorbehalten): 

'inljahres-Inhaltsverzeichnis 1931 bis 1935 der 
»chrift des Vereines Deutscher Ingenieure für 
‘5 bis 79. Format DIN A 4, 140 S. Preis 
8M. (VDI-Mitgl. 7,20 M... 


/ 


Dr.-Ing. THEODOR PÖSCHL VDI, o. Prof. a. d. 
Teehn. Hochschule Karlsruhe, Elementare 
Festiekeitslehre zum Gebrauch bei 
Vorlesuneen und zum Selbststudium 
(Lehrbuch der technischen Mechanik, zweiter Band). 
\V1/218 S. mit 156 Textabb. Berlin 1936, Verlag 
Julius Springer. Preis geb. 14,25 M. 


Dr. RÜUDOLF ROTHE, 0. Prof. a. d. Techn. Hoch- 
schule Berlin. Höhere Mathematik, Teil IV: 
Übungsaufgaben mit Lösungen, For 
melsammlung, unter Mitwirkung von >tud.- 
Rat OSKAR DEGOSANG, Assistent an der Teehn. 
Hochseh. in Berlin. 3. Heft: Integralrech- 
nung (Teubners Mathematische Leitfäden, Bi. 35). 
4) S. mit 24 Abb. Leipzig und Berlin 1936, Verlag 
B. G. Teubner. Preis kart. 2M. 


GEORGE DE BOTHEZAT, Back to Newton, 
a Challenge to Einsteins Theory of 
Relativity. V11/152 8. Verlag G. J. Stechert & 
Co. New York 1936. Preis 2,50 Dollar, 


ERNST ZIMMER, Umsturz im Weltbild 
der Physik. mit einem Geleitwort von Dr. Max 
Planck. 23. Aufl. 270 8. mit 50 Abb. München 
1936. Verlag Knorr & Hirth, G.m.b.H. Preis geb. 
5,70 


R. VON MISES, Professor an der Universität 
Istanbul, Wahrscheinliehkeit, Statistik 
und Wahrheit, Einführung in die neue Wahr- 
scheinlichkeitslehre und ihre Anwendung. Zweite 
neubearbeitete Auflage. Schriften zur wissensech. 
Weltauffassung, herausgegeben von Ph. Frank 
und M. Schlink. Bd. 3. V111/282 Wien 1936, 
Verlag ‚Julius Springer. Preis 16 M. 


Dr. PASUCUAL JQRDAN, Professor an der Univer- 
sität Rostock, Anschauliche Quanten- 
theorie. Eine Einführung in die moderne Auf- 
fassune der  Quantenerscheinungen. >. 
Berlin 1936, Verlag Julius Springer. Preis geb. 
13.80 RM. 


Dr.-Ing. ERNST SCHMIDT, Professor an der 
Technischen Hochschule Danzie, Einführung 
in die Teehnische Thermodynamik. 
V111/314 S. mit 182 Abb. und 2 Tafeln. Berlin 1936, 
Verlag Julius Springer. Preis geb. 15 M. 


Dr.-Ing. ©. FÖPPL, Professor an der Technischen 
Hochschule Braunschweig, Aufschaukelung 
und Dämpfung von Sehwingungen. 
Grundzüge der Technischen Sehwingzungslehre. 
Bd. Il. V1/212 S. mit 72 Abb. Berlin 1936, Verlag 
Julius Springer. Preis geb. 8.40 M. 


Dr. 6. HOHEISEL, Professor an der Universität 
Greifswald. Samm- 
lung Göschen Bd. 1099. 136 8. serlin und 
Leipzig 1936, Verlag W. de Gruyter & Co. Preis 
1.62 M. 


(4. MAHLER, Professor «der Mathematik und 
Physik am Gymnasium in Ulm, Physikalische 
Aufgabensammlune. Neubearbeitet von 
Prof. K. Mahler, Studienrat an der Oberreal- 
schule Aalen in Württemberg. Fünfte verbesserte 
Auflage. Sammlunz Göschen Bd. 243. 128 
Berlin und Leipzig 1936, Verlag W. de Gruyter 
& Co. Preis geb. 1,62 M. 


Dr. Mont.-Ing. F. HÖNIG, Leoben, Grund- 
sesetze der Zerkleinerung  \VDI-For- 
sechungsheft, Nr. 378. 21 8. mit 45 Abb. Berlin 1036, 
Verlag des Vereins Deutscher Ingenieure. Preis 


>M. 
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Nachriehten 


Ztsehr. f. ang 
Math. undiM: 


First Report on Viscosity and Plastieity prepared 
by the Committee for the Study of Viscosity of 
the Academy of Seiences at Amsterdam. V1/256 S. 
Amsterdam 1935 N.V.. Noord-Hollandsche  Uit- 
zevers-Maatschappij. Preis geb. 11 fl. 


H. SOLBERG, Über die freien 
euneen einer homogenen Flüssig- 
keitsschieht auf der rotierenden 
(Astrophysica Norvegiea, Vol. I, No. 7), von 
S. 237 bis 340, mit 4 Textfiguren. Oslo 1936. 


Schwin- 


Dr. GERHARD HERZBERG, V’rofessor für Physik 
an der T. H. Darmstadt, z. Zt. Carnegie Gastpro- 
fessor für Physik an der Universität Saskatsche- 
wan, Kanada, Atomspektren und Atom- 
struktur. Eine Einführung für Chemiker, Phy- 
siker und Physiko-Chemiker.  Wissenschaftliche 
Forschungsberiehte, Naturwissenschaftliche Reihe, 


herauszereben von Dr. R. E. Liesegang, Frankfurt 


a.M. Band 37. XV] ISS S. mit 70 Abb. und 21 
Tabellen. Dresden und Leipzig 1936, Verlag Th. 


Steinkopff. Preis geb. 14 M. 


Dr. S. FLÜGGE, Institut für theoretische Phy: 
an der Universität Leipzig, und Dr. A. KREE ;, 
Institut für physikalische Grundlagen der Medi: ı 
an der Universität Frankfurt a. M. Expeı 
mentelle Grundlagen der Welle 
mechanik. Wissensch. Forschungsberichte, \ 
turwissensch. Reihe, herausgegeben von Dr. R. 
Liesegang, Frankfurt a. M. Band 38. X/236 
mit 92 Abb. Dresden u. Leipzig 1936, Verlag "7. 
Steinkopff. Preis 17 


Dr. EMIL MÜLLER, Lehrbuch der dar. 
stellenden Geometrie, vollständig neu be- 
arbeitet von Dr. Erwin Kruppa, Professor an 
der Technischen Hochschule in Wien. Vierte Auf- 
lage in drei Teilen. — Teil I, Projektion auf eine 


Bildebene. VI1/130 8. mit 135 Abb. Preis kart, 
M. — Teil 1, Zugeordnete Normalrisse, 


Krumme Flächen. VI1/112 S. mit 101 Abb. Preis 
kart. 6,500 M. — Teil IH. Axonometrie, Perspektivv. 
Landkartenentwürfe. V1V148 8. mit 130 Ab). 
P’reis kart. 8,60 M. Leipzig und Berlin 1936, Verlag 
B. G. Teubner. 


NACHRICHTEN 


Prolessor Dr. Werner Friedli gestorben. 

A\m 14. September 1936 ist nach kurzer schwerer 
Krankheit im Alter von erst 43 Jahren der Pro- 
fessor für mathematisch-teehnische Versiecherungs- 
wissenschaft an der Universität Bern, Dr, 
Werner Friedli. gestorben. Kin Schüler 
Mosers hat er (die von seinem Lehrer und 
vorher von Kinkelin in der Schweiz geschaffene 
eifriee der  Versicherungsmathematik in 


hervorragendem Mabe fortgesetzt. wie seine 
seiner Sehüler Arbeiten bekunden. die zumeist 
in den Mitteilungen «der Vereinigung  Schweize- 


rischer Versicherungsmathematiker erschienen sind. 
Sie behandeln überwiegend Fragen der Pensions- 
umd Sozialversicherung. Wiederholt hat er sich auch 
über methodische Fragen der Versicherungsmathe- 
matik allgemein geäußert. wobei er anknüpfend an 
Mosersche Arbeiten (lie eleutung der Intensitäts- 
funktion hervorhebt. In einer andern Arbeit „Prä- 
zision und Approximation in der Versicherungslehre“ 
wendet er sich mit Recht gegen (die übermäßige 
Stellenzahl bei Sterbetafeln. In seiner letzten Ar- 
beit, die Anfang «des Sommers 1936 erschienen ist. 
leitet er eine einfache Beziehung zwischen den Mo- 
menten beim Gaubßschen Fehlerzesetz ab: 

3 

2.M, AM, 
Sie zeigt sich schon bei einem Material von ge- 
ringem Umfang (Zahl der Gewittertage im Jahre 1934 
bei 14 schweizerischen Stationen) recht gut be- 
stätigt. Allerdings bliebe noch, wie Friedli selbst 
sart, zu untersuchen. ob diese Beziehung für das 
Gaußsche Fehlerzesetz kennzeichnend ist. Als Prak- 
tiker war Professor Friedli im Amt für Schweize- 
rische Sozialversicherung tätie, Er war auch Mit- 
xlied des Beirats der \ereinizunge Schweizerischer 
Versicherungsmathematiker, in ausführlicher 
Nachruf wird in den Blättern für Versiecherungs- 
mathematik erscheinen. Lhorev. 


| 


Gesellschaft für angewandte Mathematik u. Mechanik. 
Hauptversammlung in Dresden. 

Im Rahmen der 94. Versammlung der Gesellschaft 

deutscher Naturforscher und Ärzte fand in der Zeit 


vom 24. bis zum 26. September in Dresden die or- 
dentliche Hauptversammlung und wissenschaftliche 
Tagung der Gesellschaft tür angewandte Mathe- 
matik und Mechanik statt. Die Vorträge des ersten 
Tages waren in der Hauptsache Fragen der Elasti- 
zitätstheorie gewidmet. Es sprachen die Herren 
Schapitz, Marguerre, Gran Olsson, Ludwig, Sör- 
rensen, Schultz-Grunow (Der Spannungsverlauf in 
rotierenden Scheiben mit exzentrischen Löchern) unil 
Weber '). Am Freitag, 25. 9. fand zunächst unter 
l.eitunze von Herrn Prandtl die ordentliche Haupt- 
versammlung der GAMM statt, über «deren Verlauf 
der ausführliche Jahresbericht die Mitglieder der 
Gesellschaft unterriehten wird. Sodann behandelten 
die Herren Prandtl. Kinner, Schliehting (Tragflügel- 
theorie bei UÜberschallgesehwindigkeit), Schmilt. 
Weinige, Pöschl, Tollmien, Weinblum und von Ohne- 
sorge ') in der Hauptsache Fragen der Strömungs- 
forschung. Außerdem führte Herr Krutzsch 
interessante kinematographische Aufnahmen von 
stabilen und instabilen Wirbelringen vor, die im 
l’hysikalischen Institut der technischen Hochschule 
Dres.ten hergestellt waren. Der letzte Tag brachte 
die Vorträge der Herren Beyer, Vogrelpohl und Föt- 
tinger, der für Herrn Henschke ') sprach. 


sehr 


Hauptversammlung wie Vorträge waren auber- 
ordentlich gut besucht. Besichtigungen von Inst'- 
tuten der technischen Hochschule schlossen sich 
den Sitzungen an.  Gesellige Zusammenkünite 
saben den Anwesenden Gelegenheit sich gegt- 
seitir kennen zu lernen, Wohl sämtliche Teilne' 
mer waren von dem Verlauf der Tagung sehr \v 
friedigt. Vor allem war man erfreut über die gro! 
Zahl allgemein interessierender Vorträge; man | 
dauerte nur, daß Herr Trefftz, der die Tagung © 
berufen und im wesentlichen vorbereitet ha! 
dureh eine schwere Erkrankung verhindert war, © 
der Tagung selbst teilzunehmen. 


Willers. 


I) Die hier nieht angeführten Vortragsthemen 
man auf Seite 256 des vorigen Heftes dieser Zeitschı 
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